<!DOCTYPE html><html lang="en"><head><meta charset="UTF-8"><meta name="viewport" content="width=device-width,initial-scale=1,maximum-scale=2"><meta name="theme-color" content="#FFF"><meta name="baidu-site-verification" content="code-DvHFUp3jdk"><meta name="baidu-site-verification" content="code-eHnRZZHAeP"><meta name="baidu-site-verification" content="code-EmlZBboqEW"><meta name="baidu-site-verification" content="code-Hf9ZXEHZ41"><meta name="360-site-verification" content="941073b91e9afaf02748ba000ba47e80"><meta name="360-site-verification" content="010bd45f0d523ce46f298dd597ef0289"><meta name="sogou_site_verification" content="hlLPNfmIjv"><meta name="sogou_site_verification" content="y3La2dTB20"><meta name="msvalidate.01" content="03829B28CC46A7AA4B5503897F475BA1"><script>!function(){var s="https://s.ssl.qhres2.com/ssl/ab77b6ea7f3fbf79.js";document.write('<script src="'+s+'" id="sozz"><\/script>')}()</script><script>!function(t,e,c,n,r,a,s){t[c]=t[c]||function(){(t[c].q=t[c].q||[]).push(arguments)},a=e.createElement(n),a.async=1,a.src="https://www.clarity.ms/tag/"+r+"?ref=bwt",s=e.getElementsByTagName(n)[0],s.parentNode.insertBefore(a,s)}(window,document,"clarity","script","8xwpfccwj2")</script><script>var _hmt=_hmt||[];!function(){var e=document.createElement("script");e.src="https://hm.baidu.com/hm.js?1d318b0509ca0909017253ed6f659ae4";var t=document.getElementsByTagName("script")[0];t.parentNode.insertBefore(e,t)}()</script><script>var _hmt=_hmt||[];!function(){var e=document.createElement("script");e.src="https://hm.baidu.com/hm.js?aea6c1471d6299bdd1c87a0836f18675";var t=document.getElementsByTagName("script")[0];t.parentNode.insertBefore(e,t)}()</script><script>var _hmt=_hmt||[];!function(){var e=document.createElement("script");e.src="https://hm.baidu.com/hm.js?498ca413655655a62058d6b5d7c62786";var t=document.getElementsByTagName("script")[0];t.parentNode.insertBefore(e,t)}()</script><script>var _hmt=_hmt||[];!function(){var e=document.createElement("script");e.src="https://hm.baidu.com/hm.js?e19e8161e13f718806b79f49240583b4";var t=document.getElementsByTagName("script")[0];t.parentNode.insertBefore(e,t)}()</script><link rel="apple-touch-icon" sizes="180x180" href="/images/apple-touch-icon.png"><link rel="icon" type="image/ico" sizes="32x32" href="/images/favicon.ico"><meta http-equiv="Cache-Control" content="no-transform"><meta http-equiv="Cache-Control" content="no-siteapp"><meta name="google-site-verification" content="UA-139249935-1"><meta name="baidu-site-verification" content="cOtBu9CtTEIuwJeN"><link rel="alternate" type="application/rss+xml" title="了尘兰若的小坑" href="https://liaochenlanruo.gitee.io/rss.xml"><link rel="alternate" type="application/atom+xml" title="了尘兰若的小坑" href="https://liaochenlanruo.gitee.io/atom.xml"><link rel="alternate" type="application/json" title="了尘兰若的小坑" href="https://liaochenlanruo.gitee.io/feed.json"><link rel="stylesheet" href="//fonts.googleapis.com/css?family=Mulish:300,300italic,400,400italic,700,700italic%7CFredericka%20the%20Great:300,300italic,400,400italic,700,700italic%7CNoto%20Serif%20JP:300,300italic,400,400italic,700,700italic%7CNoto%20Serif%20SC:300,300italic,400,400italic,700,700italic%7CInconsolata:300,300italic,400,400italic,700,700italic&display=swap&subset=latin,latin-ext"><link rel="stylesheet" href="/css/app.css?v=0.2.5"><meta name="keywords" content="比较基因组学"><link rel="canonical" href="https://liaochenlanruo.gitee.io/post/848f.html"><title>PGCGAP中文说明 - 生物信息 | ResearchGo = 了尘兰若的小坑 = liaochenlanruo</title><meta name="generator" content="Hexo 5.4.0"></head><body itemscope itemtype="http://schema.org/WebPage"><div id="loading"><div class="cat"><div class="body"></div><div class="head"><div class="face"></div></div><div class="foot"><div class="tummy-end"></div><div class="bottom"></div><div class="legs left"></div><div class="legs right"></div></div><div class="paw"><div class="hands left"></div><div class="hands right"></div></div></div></div><div id="container"><header id="header" itemscope itemtype="http://schema.org/WPHeader"><div class="inner"><div id="brand"><div class="pjax"><h1 itemprop="name headline">PGCGAP中文说明</h1><div class="meta"><span class="item" title="Created: 2019-04-28 15:22:33"><span class="icon"><i class="ic i-calendar"></i> </span><span class="text">Posted on</span> <time itemprop="dateCreated datePublished" datetime="2019-04-28T15:22:33+08:00">2019-04-28</time> </span><span class="item" title="Symbols count in article"><span class="icon"><i class="ic i-pen"></i> </span><span class="text">Symbols count in article</span> <span>33k</span> <span class="text">words</span> </span><span class="item" title="Reading time"><span class="icon"><i class="ic i-clock"></i> </span><span class="text">Reading time</span> <span>30 mins.</span></span></div></div></div><nav id="nav"><div class="inner"><div class="toggle"><div class="lines" aria-label="Toggle navigation bar"><span class="line"></span> <span class="line"></span> <span class="line"></span></div></div><ul class="menu"><li class="item title"><a href="/" rel="start">ResearchGo</a></li></ul><ul class="right"><li class="item theme"><i class="ic i-sun"></i></li><li class="item search"><i class="ic i-search"></i></li></ul></div></nav></div><div id="imgs" class="pjax"><img src="https://fastly.jsdelivr.net/gh/liaochenlanruo/cdn@master/img/custom/bgs/thumb_81.webp"></div></header><div id="waves"><svg class="waves" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink" viewBox="0 24 150 28" preserveAspectRatio="none" shape-rendering="auto"><defs><path id="gentle-wave" d="M-160 44c30 0 58-18 88-18s 58 18 88 18 58-18 88-18 58 18 88 18 v44h-352z"/></defs><g class="parallax"><use xlink:href="#gentle-wave" x="48" y="0"/><use xlink:href="#gentle-wave" x="48" y="3"/><use xlink:href="#gentle-wave" x="48" y="5"/><use xlink:href="#gentle-wave" x="48" y="7"/></g></svg></div><main><div class="inner"><div id="main" class="pjax"><div class="article wrap"><div class="breadcrumb" itemscope itemtype="https://schema.org/BreadcrumbList"><i class="ic i-home"></i> <span><a href="/">Home</a></span><i class="ic i-angle-right"></i> <span class="current" itemprop="itemListElement" itemscope itemtype="https://schema.org/ListItem"><a href="/categories/Bioinformatics/" itemprop="item" rel="index" title="In 生物信息"><span itemprop="name">生物信息</span></a><meta itemprop="position" content="1"></span></div><article itemscope itemtype="http://schema.org/Article" class="post block" lang="en"><link itemprop="mainEntityOfPage" href="https://liaochenlanruo.gitee.io/post/848f.html"><span hidden itemprop="author" itemscope itemtype="http://schema.org/Person"><meta itemprop="image" content="/images/head.jpg"><meta itemprop="name" content="Hualin Liu"><meta itemprop="description" content="liaochenlanruo, 分享微生物生物信息学分析方法，欢迎加入QQ群交流945751012，不接受群内广告！"></span><span hidden itemprop="publisher" itemscope itemtype="http://schema.org/Organization"><meta itemprop="name" content="了尘兰若的小坑"></span><div class="body md" itemprop="articleBody"><div class="gallery" itemscope itemtype="http://schema.org/ImageGallery"><img data-src="https://fastly.jsdelivr.net/gh/liaochenlanruo/cdn@master/img/custom/bgs/thumb_81.webp" itemprop="contentUrl"></div><p><strong>为了方便广大中文用户学习 PGCGAP 的使用，特意书写该中文文档，但中文文档更新较慢，强烈建议大家阅读英文文档！</strong></p><span id="more"></span><p><span class="exturl" data-url="aHR0cHM6Ly9saWFvY2hlbmxhbnJ1by5naXRodWIuaW8vcGdjZ2Fw">English Readme</span> | <span class="exturl" data-url="aHR0cHM6Ly9saWFvY2hlbmxhbnJ1by5naXRodWIuaW8vMjAxOS8wNC8yOC9QR0NHQVAlRTQlQjglQUQlRTYlOTYlODclRTglQUYlQjQlRTYlOTglOEUv">中文说明</span></p><h2 id="简介"><a class="anchor" href="#简介">#</a> 简介</h2><p>PGCGAP 是用于原核生物基因组学和比较基因组学分析管道，目前该管道包含 12 个模块，可以接受 Illumina 双端 reads、Oxford reads 或 PacBio reads 作为输入，可以完成基因组组装、基因预测和注释，并可以进行比较基因组学分析，包括构建单拷贝核心蛋白进化树以及单拷贝核心基因 SNPs 进化树，泛基因组分析与进化树构建，全基因组平均核苷酸一致性（ANI）计算，同源蛋白家族聚类及进化树构建，COG 注释，SNPs 和 INDELs calling，抗生素抗性基因 / 毒力因子预测，Multi-FASTA 进化树构建，组装后基因组短序列过滤与统计信息呈现（genome size，GC content……）。</p><h2 id="安装"><a class="anchor" href="#安装">#</a> 安装</h2><p>PGCGAP 可以安装于 Windows 子系统 Linux（WSL）、Linux x64 系统以及 macOS 中。</p><p><strong>Step1：通过<span class="exturl" data-url="aHR0cHM6Ly9iaW9jb25kYS5naXRodWIuaW8vdXNlci9pbnN0YWxsLmh0bWwjaW5zdGFsbC1jb25kYQ=="> Bioconda</span> 安装 PGCGAP</strong></p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$conda</span> create -n pgcgap <span class="token assign-left variable">python</span><span class="token operator">=</span><span class="token number">3.7</span></pre></td></tr><tr><td data-num="2"></td><td><pre></pre></td></tr><tr><td data-num="3"></td><td><pre><span class="token variable">$conda</span> activate pgcgap</pre></td></tr><tr><td data-num="4"></td><td><pre></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$conda</span> <span class="token function">install</span> pgcgap</pre></td></tr><tr><td data-num="6"></td><td><pre></pre></td></tr><tr><td data-num="7"></td><td><pre><span class="token variable">$conda</span> deactivate</pre></td></tr></table></figure><p><strong>声明：用conda安装时一直在“Solving environment”怎么办？</strong>随着加入conda的软件越来越多，使其索引库变的庞大，因此安装新的软件时需要逐一验证软件间的兼容性，会耗费大量时间。严重的时候会无法完成软件的安装，这个时候坐以待毙是不行的。下面介绍两种办法解决此问题：<br></p><ul><li><p>Method 1: 使用 mamba (非常快) 替代 conda。前提是已经用 conda 创建好了 pgcgap 的虚拟环境：</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$conda</span> activate pgcgap</pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$conda</span> <span class="token function">install</span> mamba -c conda-forge</pre></td></tr><tr><td data-num="3"></td><td><pre><span class="token variable">$mamba</span> <span class="token function">install</span> pgcgap</pre></td></tr></table></figure></li><li><p>Method 2: 使用本人提供的 pgcgap<span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL2xpYW9jaGVubGFucnVvL3BnY2dhcC9ibG9iL21hc3Rlci9jb25kYS9wZ2NnYXBfbGF0ZXN0X2Vudi55bQ=="> 配置文件</span>来创建环境并安装 PGCGAP：</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment"># download pgcgap_latest_env.yml</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$wget</span> https://github.com/liaochenlanruo/pgcgap/blob/master/conda/pgcgap_latest_env.yml</pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment"># create a conda environment named as pgcgap and install the latest version of PGCGAP</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$conda</span> <span class="token function">env</span> create -f pgcgap_latest_env.yml</pre></td></tr></table></figure></li></ul><br><p><strong>Step2：配置 COG 数据库</strong> (初次安装 PGCGAP 后需要执行此步骤)</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$conda</span> activate pgcgap</pre></td></tr><tr><td data-num="2"></td><td><pre></pre></td></tr><tr><td data-num="3"></td><td><pre><span class="token variable">$pgcgap</span> --setup-COGdb</pre></td></tr><tr><td data-num="4"></td><td><pre></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$conda</span> deactivate</pre></td></tr></table></figure><p><strong>Step3: 升级 PGCGAP (升级版本时运行)</strong></p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$conda</span> activate pgcgap</pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$conda</span> update pgcgap</pre></td></tr><tr><td data-num="3"></td><td><pre><span class="token comment"># v1.0.28 以后可通过如下命令升级</span></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token variable">$pgcgap</span> --check-update</pre></td></tr></table></figure><p><strong>此外，用户也可以通过容器 (docker) 安装 PGCGAP</strong></p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$docker</span> pull quay.io/biocontainers/pgcgap:<span class="token operator">&lt;</span>tag<span class="token operator">></span></pre></td></tr></table></figure><p><strong>注：</strong> 前提是用户电脑中安装了<span class="exturl" data-url="aHR0cHM6Ly93d3cuZG9ja2VyLmNvbS8="> Docker</span>，Docker 可以跨平台使用。可用的 tag 可<span class="exturl" data-url="aHR0cHM6Ly9xdWF5LmlvL3JlcG9zaXRvcnkvYmlvY29udGFpbmVycy9wZ2NnYXA/dGFiPXRhZ3M=">在此</span>查询，建议安装最新版。</p><h2 id="依赖包"><a class="anchor" href="#依赖包">#</a> 依赖包</h2><ul><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL2FicmljYXRl">Abricate</span></li><li><span class="exturl" data-url="aHR0cDovL3d3dy5iY2dzYy5jYS9wbGF0Zm9ybS9iaW9pbmZvL3NvZnR3YXJlL2FieXNzLw==">ABySS</span></li><li><span class="exturl" data-url="aHR0cDovL2NhbnUucmVhZHRoZWRvY3Mub3JnLw==">Canu</span></li><li><span class="exturl" data-url="aHR0cDovL3dlaXpob25nbGktbGFiLm9yZy9jZC1oaXQv">CD-HIT</span></li><li><span class="exturl" data-url="aHR0cHM6Ly93d3cuZ251Lm9yZy9zb2Z0d2FyZS9jb3JldXRpbHMv">Coreutils</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL2JidWNoZmluay9kaWFtb25k">Diamond</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL1BhckJMaVNTL0Zhc3RBTkk=">FastANI</span></li><li><span class="exturl" data-url="aHR0cDovL3d3dy5hdGdjLW1vbnRwZWxsaWVyLmZyL2Zhc3RtZS9iaW5hcmllcy5waHA=">Fastme</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL09wZW5HZW5lL2Zhc3Rw">Fastp</span></li><li><span class="exturl" data-url="aHR0cDovL3d3dy5taWNyb2Jlc29ubGluZS5vcmcvZmFzdHRyZWU=">FastTree</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3Nhbmdlci1wYXRob2dlbnMvZ3ViYmlucw==">Gubbins</span> &gt;=2.3.4</li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3NhbXRvb2xzL2h0c2xpYg==">Htslib</span></li><li><span class="exturl" data-url="aHR0cDovL3d3dy5pcXRyZWUub3JnLw==">IQ-TREE</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9tYWZmdC5jYnJjLmpwL2FsaWdubWVudC9zb2Z0d2FyZS8=">Mafft</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL21hcmJsL01hc2g=">Mash</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3NvZWRpbmdsYWIvbW1zZXFzMg==">Mmseqs2</span></li><li><span class="exturl" data-url="aHR0cHM6Ly93d3cuZWJpLmFjLnVrL1Rvb2xzL21zYS9tdXNjbGUv">Muscle</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9ibGFzdC5uY2JpLm5sbS5uaWguZ292L0JsYXN0LmNnaT9DTUQ9V2ViJmFtcDtQQUdFX1RZUEU9Qmxhc3REb2NzJmFtcDtET0NfVFlQRT1Eb3dubG9hZA==">NCBI-blast+</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL2RhdmlkZW1tcy9PcnRob0ZpbmRlcg==">OrthoFinder</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9vcGVuamRrLmphdmEubmV0Lw==">OpenJDK8</span></li><li><span class="exturl" data-url="aHR0cDovL3d3dy5ib3JrLmVtYmwuZGUvcGFsMm5hbC8=">PAL2NAL v14</span></li><li><span class="exturl" data-url="aHR0cDovL3RyaW1hbC5jZ2Vub21pY3Mub3Jn">trimAL</span></li><li><span class="exturl" data-url="aHR0cDovL3d3dy5wZXJsLm9yZy9nZXQuaHRtbA==">Perl</span> &amp; the modules<ul><li><span class="exturl" data-url="aHR0cDovL21ldGFjcGFuLm9yZy9wb2QvQmlvUGVybA==">perl-bioperl</span></li><li><span class="exturl" data-url="aHR0cDovL21ldGFjcGFuLm9yZy9wb2QvRGF0YTo6RHVtcGVy">perl-data-dumper</span></li><li><span class="exturl" data-url="aHR0cDovL21ldGFjcGFuLm9yZy9wb2QvRmlsZTo6VGVl">perl-file-tee</span></li><li><span class="exturl" data-url="aHR0cDovL21ldGFjcGFuLm9yZy9wb2QvR2V0b3B0OjpMb25n">perl-getopt-long</span></li><li><span class="exturl" data-url="aHR0cDovL3NlYXJjaC5jcGFuLm9yZy9+bWFyZWtyL1BvZC1Vc2FnZS0xLjY5Lw==">perl-pod-usage</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9tZXRhY3Bhbi5vcmcvcG9kL3JlbGVhc2UvRExVWC9QYXJhbGxlbC1Gb3JrTWFuYWdlci0wLjcuNS9Gb3JrTWFuYWdlci5wbQ==">perl-parallel-forkmanager</span></li></ul></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL3Byb2trYQ==">Prokka</span></li><li><span class="exturl" data-url="aHR0cHM6Ly93d3cucHl0aG9uLm9yZy8=">Python</span> &amp; the modules<ul><li><span class="exturl" data-url="aHR0cHM6Ly9iaW9weXRob24ub3JnLw==">biopython</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9tYXRwbG90bGliLm9yZy8=">matplotlib</span></li><li><span class="exturl" data-url="aHR0cDovL3d3dy5udW1weS5vcmcv">numpy</span></li><li><span class="exturl" data-url="aHR0cDovL3BhbmRhcy5weWRhdGEub3JnLw==">pandas</span></li><li><span class="exturl" data-url="aHR0cDovL3NlYWJvcm4ucHlkYXRhLm9yZy8=">seaborn</span></li></ul></li><li><span class="exturl" data-url="aHR0cHM6Ly93d3cuci1wcm9qZWN0Lm9yZy8=">R</span> &amp; the packages<ul><li><span class="exturl" data-url="aHR0cHM6Ly9jcmFuLnItcHJvamVjdC5vcmcvd2ViL3BhY2thZ2VzL2NvcnJwbG90L2luZGV4Lmh0bWw=">corrplot</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9jcmFuLnItcHJvamVjdC5vcmcvd2ViL3BhY2thZ2VzL2dncGxvdDIv">ggplot2</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9jcmFuLnItcHJvamVjdC5vcmcvd2ViL3BhY2thZ2VzL2dwbG90cy8=">gplots</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9jcmFuLnItcHJvamVjdC5vcmcvd2ViL3BhY2thZ2VzL3BoZWF0bWFwL2luZGV4Lmh0bWw=">pheatmap</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9jcmFuLnItcHJvamVjdC5vcmcvd2ViL3BhY2thZ2VzL3Bsb3RyaXgv">plotrix</span></li></ul></li><li><span class="exturl" data-url="aHR0cHM6Ly9zYW5nZXItcGF0aG9nZW5zLmdpdGh1Yi5pby9Sb2FyeS8=">Roary</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL25ham9zaGkvc2lja2xl">Sickle-trim</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL3NuaXBweQ==">Snippy</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3Nhbmdlci1wYXRob2dlbnMvc25wLXNpdGVz">Snp-sites</span></li><li><span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3Jyd2ljay9VbmljeWNsZXI=">unicycler</span></li><li><span class="exturl" data-url="aHR0cHM6Ly93d3cuZ251Lm9yZy9zb2Z0d2FyZS93Z2V0Lw==">wget</span></li></ul><h2 id="pgcgap用法"><a class="anchor" href="#pgcgap用法">#</a> PGCGAP 用法</h2><ul><li><p><strong>显示帮助信息：</strong></p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --help</pre></td></tr></table></figure></li><li><p><strong>管道用法：</strong></p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> <span class="token punctuation">[</span>modules<span class="token punctuation">]</span> <span class="token punctuation">[</span>options<span class="token punctuation">]</span></pre></td></tr></table></figure></li><li><p><strong>展示各模块的参数：</strong></p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> <span class="token punctuation">[</span>Assemble<span class="token operator">|</span>Annotate<span class="token operator">|</span>ANI<span class="token operator">|</span>AntiRes<span class="token operator">|</span>CoreTree<span class="token operator">|</span>MASH<span class="token operator">|</span>OrthoF<span class="token operator">|</span>Pan<span class="token operator">|</span>pCOG<span class="token operator">|</span>VAR<span class="token operator">|</span>STREE<span class="token operator">|</span>ACC<span class="token punctuation">]</span></pre></td></tr></table></figure></li><li><p><strong>展示各模块的运行示例：</strong>（这货是我用的最多的）</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> Examples</pre></td></tr></table></figure></li><li><p><strong>配置 COG 数据库：</strong> (初次安装 PGCGAP 后需要配置 COG 数据库)</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --setup-COGdb</pre></td></tr></table></figure></li><li><p><strong>功能模块：</strong></p><ul><li><p><strong>[--All]</strong> 运行 Assemble, Annotate, CoreTree, Pan, OrthoF, ANI, MASH 和 pCOG 模块</p></li><li><p><strong>[--Assemble]</strong> 基因组组装</p></li><li><p><strong>[--Annotate]</strong> 基因预测及注释</p></li><li><p><strong>[--CoreTree]</strong> 构建单拷贝核心蛋白进化树与核心 SNPs 进化树</p></li><li><p><strong>[--Pan]</strong> 泛基因组分析并构建单拷贝核心蛋白进化树</p></li><li><p><strong>[--OrthoF]</strong> 同源蛋白家族聚类及单拷贝直系同源蛋白进化树构建</p></li><li><p><strong>[--ANI]</strong> 计算平均核苷酸一致性 (ANI)</p></li><li><p><strong>[--MASH]</strong> 通过 MinHash 估算基因组 / 宏基因组相似性</p></li><li><p><strong>[--pCOG]</strong> COG 注释</p></li><li><p><strong>[--VAR]</strong> 变异检测并构建核心基因组进化树</p></li><li><p><strong>[--AntiRes]</strong> 从基因组 (contigs/scaffolds) 中预测抗生素抗性基因或毒力基因</p></li><li><p><strong>[--STREE]</strong> 基于 Multi-FASTA 序列 (所有序列在一个文件中) 构建系统发育树</p></li><li><p><strong>[--ACC]</strong> 一些实用的附加程序 (目前只开发了 &quot;Assess&quot; 用于对基因组中的短序列进行过滤，并评估过滤前后的基因组状态)</p></li></ul></li><li><p><strong>全局参数（请参照英文版，参数有所改变，中文版暂时没有时间修改）：</strong></p><ul><li><p><strong>[--strain_num (INT)]</strong> [Required by &quot;--All&quot;, &quot;--CoreTree&quot;, &quot;--Pan&quot;, &quot;--VAR&quot; and &quot;--COG&quot;] 用于分析的菌株数目，不包含参考基因组</p></li><li><p><strong>[--ReadsPath (PATH)]</strong> [Required by &quot;--All&quot;, &quot;--Assemble&quot; and &quot;--VAR&quot;] 所有菌株测序 reads 所在的目录路径 Default ./Reads/Illumina)</p></li><li><p><strong>[--scafPath (PATH)]</strong> [Required by &quot;--All&quot;, &quot;--Assess&quot;, &quot;--Annotate&quot; and &quot;--MASH&quot;] contigs/scaffolds 的存放路径 (Default &quot;Results/Assembles/Scaf/Illumina&quot;)</p></li><li><p><strong>[--AAsPath (PATH)]</strong> [Required by &quot;--All&quot;, &quot;--CoreTree&quot;, &quot;--OrthoF&quot; and &quot;--pCOG&quot;] 所有菌株的氨基酸序列文件的存放路径 (Default &quot;./Results/Annotations/AAs&quot;)</p></li><li><p><strong>[--reads1 (STRING)]</strong> [Required by &quot;--All&quot;, &quot;--Assemble&quot; and &quot;--VAR&quot;] reads 1 的后缀名 (例如 reads 1 的名字为 &quot;YBT-1520_L1_I050.R1.clean.fastq.gz&quot;，&quot;YBT-1520&quot; 是菌株名，则后缀名为 &quot;.R1.clean.fastq.gz&quot;)</p></li><li><p><strong>[--reads2 (STRING)]</strong> [Required by &quot;--All&quot;, &quot;--Assemble&quot; and &quot;--VAR&quot;] reads 2 的后缀名</p></li><li><p><strong>[--Scaf_suffix (STRING)]</strong> [Required by &quot;--All&quot;, &quot;--Assess&quot;, &quot;--Annotate&quot; &quot;MASH&quot; and &quot;--ANI&quot;] contigs/scaffolds 的后缀名 (Default -8.fa)</p></li><li><p><strong>[--filter_length (INT)]</strong> [Required by &quot;--All&quot;, &quot;--Assemble&quot; and &quot;--Assess&quot;]&gt; Sequences shorter than the 'filter_length' will be deleted from the assembled genomes. ( Default 200 )</p></li><li><p><strong>[--codon (INT)]</strong> [Required by &quot;--All&quot;, &quot;--Annotate&quot;, &quot;--CoreTree&quot; and &quot;--Pan&quot;] 翻译密码子表 (Default 11)</p><ul><li>1 Universal code</li><li>2 Vertebrate mitochondrial code</li><li>3 Yeast mitochondrial code</li><li>4 Mold, Protozoan, and Coelenterate Mitochondrial code and Mycoplasma/Spiroplasma code</li><li>5 Invertebrate mitochondrial</li><li>6 Ciliate, Dasycladacean and Hexamita nuclear code</li><li>9 Echinoderm and Flatworm mitochondrial code</li><li>10 Euplotid nuclear code</li><li>11 Bacterial, archaeal and plant plastid code ( Default )</li><li>12 Alternative yeast nuclear code</li><li>13 Ascidian mitochondrial code</li><li>14 Alternative flatworm mitochondrial code</li><li>15 Blepharisma nuclear code</li><li>16 Chlorophycean mitochondrial code</li><li>21 Trematode mitochondrial code</li><li>22 Scenedesmus obliquus mitochondrial code</li><li>23 Thraustochytrium mitochondrial code</li></ul></li><li><p><strong>[--suffix_len (INT)]</strong> [Required by &quot;--All&quot;, &quot;--Assemble&quot; and &quot;--VAR&quot;] <strong>(强烈建议设置此项)</strong> reads 后缀的长度。例如 &quot;YBT-1520_L1_I050.R1.clean.fastq.gz&quot; 的 --suffix_len 为 26 (&quot;YBT-1520&quot; 为菌株名) (Default 0)</p></li><li><p><strong>[--logs (STRING)]</strong> Log 文件的名字 (Default Logs.txt)</p></li><li><p><strong>[--threads (INT)]</strong> 运行程序时调用的线程数目 (Default 4)<br><br></p></li></ul></li><li><p><strong>各模块的局部参数：</strong></p><ul><li><p><strong>--Assemble</strong></p><ul><li><p><strong>[--platform (STRING)]</strong> [Required] 测序平台，可以选择 “illumina”, “pacbio” 和 “oxford” (Default illumina)</p></li><li><p><strong>[--assembler (STRING)]</strong> [Required] 用于 illumina 数据组装的软件，可选 &quot;abyss&quot;, &quot;spades&quot; 或 &quot;auto&quot; ( Default abyss )</p></li><li><p><strong>[--kmmer (INT)]</strong> [Required] Illumina 数据组装时采用的 k-mer 大小 (Default 81)</p></li><li><p><strong>[--genomeSize (FLOAT)]</strong> [Required] 预估的基因组大小，如 3.7m、2.8g，组装 PacBio data 和 Oxford data 时需要设置此项 (Default Unset)</p></li><li><p><strong>[--short1 (STRING)]</strong> [Required] FASTQ file of first short reads in each pair. Needed by hybrid assembly ( Default Unset )</p></li><li><p><strong>[--short2 (STRING)]</strong> [Required] FASTQ file of second short reads in each pair. Needed by hybrid assembly ( Default Unset )</p></li><li><p><strong>[--long (STRING)]</strong> [Required] FASTQ or FASTA file of long reads. Needed by hybrid assembly ( Default Unset )</p></li><li><p><strong>[--hout (STRING)]</strong> [Required] Output directory for hybrid assembly ( Default ../../Results/Assembles/Hybrid )</p></li></ul></li><li><p><strong>--Annotate</strong></p><ul><li><p><strong>[--genus (STRING)]</strong> 菌株的属名 (Default &quot;NA&quot;)</p></li><li><p><strong>[--species (STRING)]</strong> 菌株的种名 (Default &quot;NA&quot;)<br><br></p></li></ul></li><li><p><strong>--CoreTree</strong></p><ul><li><p><strong>[--CDsPath (PATH)]</strong> [Required] 包含所有菌株核苷酸序列文件的路径，如果设置为 &quot;NO&quot;，将不会构建核心 SNPs 进化树 ( Default &quot;./Results/Annotations/CDs&quot; )</p></li><li><p><strong>[-c (FLOAT)]</strong> 序列一致性 (identity) 阈值 ( Default 0.5)</p></li><li><p><strong>[-n (INT)]</strong> Word_length, -n 2 for thresholds 0.4-0.5, -n 3 for thresholds 0.5-0.6, -n 4 for thresholds 0.6-0.7, -n 5 for thresholds 0.7-1.0 ( Default 2 )</p></li><li><p><strong>[-G (INT)]</strong> Use global (set to 1) or local (set to 0) sequence identity, ( Default 0 )</p></li><li><p><strong>[-t (INT)]</strong> Tolerance for redundance ( Default 0 )</p></li><li><p><strong>[-aL (FLOAT)]</strong> Alignment coverage for the longer sequence. If set to 0.9, the alignment must covers 90% of the sequence ( Default 0.5 )</p></li><li><p><strong>[-aS (FLOAT)]</strong> Alignment coverage for the shorter sequence. If set to 0.9, the alignment must covers 90% of the sequence ( Default 0.7 )</p></li><li><p><strong>[-g (INT)]</strong> If set to 0, a sequence is clustered to the first cluster that meet the threshold (fast cluster). If set to 1, the program will cluster it into the most similar cluster that meet the threshold (accurate but slow mode, Default 1)</p></li><li><p><strong>[-d (INT)]</strong> length of description in .clstr file. if set to 0, it takes the fasta defline and stops at first space ( Default 0 )<br><br></p></li></ul></li><li><p><strong>--Pan</strong></p><ul><li><p><strong>[--GffPath (PATH)]</strong> [Required] 存放所有菌株 GFF3 格式文件的路径 ( Default &quot;./Results/Annotations/GFF&quot; )</p></li><li><p><strong>[--identi (INT)]</strong> Minimum percentage identity for blastp ( Default 95 )<br><br></p></li></ul></li><li><p><strong>--OrthoF</strong></p><ul><li><strong>[--Sprogram (STRING)]</strong> 序列对比程序，Options: blast, mmseqs, blast_gz, diamond (Default blast)<br><br></li></ul></li><li><p><strong>--ANI</strong></p><ul><li><p><strong>[--queryL (FILE)]</strong> [Required] The file containing paths to query genomes, one per line ( Default scaf.list )</p></li><li><p><strong>[--refL (FILE)]</strong> [Required] The file containing paths to reference genomes, one per line. ( Default scaf.list )</p></li><li><p><strong>[--ANIO (FILE)]</strong> The name of output file ( Default &quot;Results/ANI/ANIs&quot; )</p></li></ul></li></ul><br><ul><li><p><strong>--VAR</strong></p><ul><li><p><strong>[--refgbk (FILE)]</strong> [Required] The full path and name of reference genome in GENBANK format ( recommended ), fasta format is also OK. For example: &quot;/mnt/g/test/ref.gbk&quot;</p></li><li><p><strong>[--qualtype (STRING)]</strong> [Required] Type of quality values (solexa (CASAVA &lt; 1.3), illumina (CASAVA 1.3 to 1.7), sanger (which is CASAVA &gt;= 1.8)). ( Default sanger )</p></li><li><p><strong>[--qual (INT)]</strong> Threshold for trimming based on average quality in a window. ( Default 20 )</p></li><li><p><strong>[--length (INT)]</strong> Threshold to keep a read based on length after trimming. ( Default 20 )</p></li><li><p><strong>[--mincov (INT)]</strong> The minimum number of reads covering a site to be considered ( Default 10 )</p></li><li><p><strong>[--minfrac (FLOAT)]</strong> The minimum proportion of those reads which must differ from the reference ( Default 0.9 )</p></li><li><p><strong>[--minqual (INT)]</strong> The minimum VCF variant call &quot;quality&quot; ( Default 100 )</p></li><li><p><strong>[--ram (INT)]</strong> Try and keep RAM under this many GB ( Default 8 )</p></li><li><p><strong>[--tree_builder (STRING)]</strong> Application to use for tree building [raxml|fasttree|hybrid] ( Default fasttree)</p></li><li><p><strong>[--iterations (INT)]</strong> Maximum No. of iterations for gubbins ( Default 5 )<br><br></p></li></ul></li><li><p><strong>--AntiRes</strong></p><ul><li><p><strong>[--db (STRING)]</strong> [Required] 用于分析的数据库，options: <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjQxNDU1MzI=">argannot</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjc3ODk3MDU=">card</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wbWMvYXJ0aWNsZXMvUE1DNTM0MzEzNi8=">ecoh</span>, ecoli_vf, <span class="exturl" data-url="aHR0cHM6Ly93d3cuYmlvcnhpdi5vcmcvY29udGVudC8xMC4xMTAxLzU1MDcwN3Yx">ncbi</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjQ3NzcwOTI=">plasmidfinder</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjI3ODI0ODc=">resfinder</span> and <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjY1Nzg1NTk=">vfdb</span>. ( Default ncbi )</p></li><li><p><strong>[--identity (INT)]</strong> [Required] Minimum %identity to keep the result, should be a number between 1 to 100. ( Default 75 )</p></li><li><p><strong>[--coverage (INT)]</strong> [Required] Minimum %coverage to keep the result, should be a number between 0 to 100. ( Default 50 )</p></li></ul></li><li><p><strong>--STREE</strong></p><ul><li><p><strong>[--seqfile (STRING)]</strong> [Required] Path of the sequence file for analysis.</p></li><li><p><strong>[--seqtype (INT)]</strong> [Required] Type Of Sequence (p, d, c for Protein, DNA, Codons, respectively). ( Default p )</p></li><li><p><strong>[--bsnum (INT)]</strong> [Required] Times for bootstrap. ( Default 1000 )</p></li></ul></li><li><p><strong>--ACC</strong></p><ul><li><strong>[--Assess (STRING)]</strong> Filter short sequences in the genome and assess the status of the genome (详细参数通过 &quot;pgcgap ACC&quot; 查看)</li></ul></li></ul></li><li><p><strong>依赖软件安装目录</strong></p><p>Not needed if they were in the environment variables path. Users can check with the &quot;--check-external-programs&quot; option for the essential programs.<br><br></p></li><li><p><strong>[--abricate-bin (PATH)]</strong> Path to abyss binary file. Default tries if abyss is in PATH;</p></li><li><p><strong>[--abyss-bin (PATH)]</strong> Path to abyss binary file. Default tries if abyss is in PATH;</p></li><li><p><strong>[--canu-bin (PATH)]</strong> Path to canu binary file. Default tries if canu is in PATH;</p></li><li><p><strong>[--cd-hit-bin (PATH)]</strong> Path to cd-hit binary file. Default tries if cd-hit is in PATH;</p></li><li><p><strong>[--fastANI-bin (PATH)]</strong> Path to the fastANI binary file. Default tries if fastANI is in PATH;</p></li><li><p><strong>[--Gblocks-bin (PATH)]</strong> Path to the Gblocks binary file. Default tries if Gblocks is in PATH;</p></li><li><p><strong>[--gubbins-bin (PATH)]</strong> Path to the run_gubbins.py binary file. Default tries if run_gubbins.py is in PATH;</p></li><li><p><strong>[--iqtree-bin (PATH)]</strong> Path to the iqtree binary file. Default tries if iqtree is in PATH;</p></li><li><p><strong>[--mafft-bin (PATH)]</strong> Path to mafft binary file. Default tries if mafft is in PATH;</p></li><li><p><strong>[--mash-bin (PATH)]</strong> Path to the mash binary file. Default tries if mash is in PATH.</p></li><li><p><strong>[--modeltest-ng-bin (PATH)]</strong> Path to the modeltest-ng binary file. Default tries if modeltest-ng is in PATH.</p></li><li><p><strong>[--muscle-bin (PATH)]</strong> Path to the muscle binary file. Default tries if muscle is in PATH.</p></li><li><p><strong>[--orthofinder-bin (PATH)]</strong> Path to the orthofinder binary file. Default tries if orthofinder is in PATH;</p></li><li><p><strong>[--pal2nal-bin (PATH)]</strong> Path to the <span class="exturl" data-url="aHR0cDovL3BhbDJuYWwucGw=">pal2nal.pl</span> binary file. Default tries if <span class="exturl" data-url="aHR0cDovL3BhbDJuYWwucGw=">pal2nal.pl</span> is in PATH;</p></li><li><p><strong>[--prodigal-bin (PATH)]</strong> Path to prodigal binary file. Default tries if prodigal is in PATH;</p></li><li><p><strong>[--prokka-bin (PATH)]</strong> Path to prokka binary file. Default tries if prokka is in PATH;</p></li><li><p><strong>[--raxml-ng-bin (PATH)]</strong> Path to the raxml-ng binary file. Default tries if raxml-ng is in PATH;</p></li><li><p><strong>[--roary-bin (PATH)]</strong> Path to the roary binary file. Default tries if roary is in PATH;</p></li><li><p><strong>[--sickle-bin (PATH)]</strong> Path to the sickle-trim binary file. Default tries if sickle is in PATH.</p></li><li><p><strong>[--snippy-bin (PATH)]</strong> Path to the snippy binary file. Default tries if snippy is in PATH;</p></li><li><p><strong>[--snp-sites-bin (PATH)]</strong> Path to the snp-sites binary file. Default tries if snp-sites is in PATH;</p></li><li><p><strong>[--unicycler-bin (PATH)]</strong> Path to the unicycler binary file. Default tries if unicycler is in PATH;</p></li></ul><br><ul><li><p><strong>配置 COG 数据库</strong></p><ul><li><strong>[--setup-COGdb]</strong> 首次安装 PGCGAP 后需要执行此步<br><br></li></ul></li><li><p>检查依赖软件包是否安装 (<strong>强烈建议在安装完 PGCGAP 之后运行此步</strong>):</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --check-external-programs</pre></td></tr></table></figure></li></ul><br><h2 id="示例"><a class="anchor" href="#示例">#</a> 示例</h2><ul><li><p><strong>Example 1:</strong> 执行所有模块，以 <em>Escherichia coli</em> 的 6 个 Illumina 双端 reads 为数据集。</p><p><strong>注</strong>：为了提高灵活性，&quot;VAR&quot; 模块需要额外添加。<br></p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --All --platform illumina --filter_length <span class="token number">200</span> --ReadsPath Reads/Illumina --reads1 _1.fastq.gz --reads2 _2.fastq.gz --suffix_len <span class="token number">11</span> --kmmer <span class="token number">81</span> --genus Escherichia --species “Escherichia coli” --codon <span class="token number">11</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --VAR --refgbk /mnt/h/PGCGAP_Examples/Reads/MG1655.gbff --qualtype sanger</pre></td></tr></table></figure></li><li><p><strong>Example 2:</strong> 基因组组装。</p><ul><li><p><strong>Illumina 双端 reads 组装</strong></p><p>该数据集中，reads 的命名格式为 “strain_1.fastq.gz” 和 “strain_2.fastq.gz”。 后缀名为 “_1.fastq.gz”，其长度为 11，因此 &quot;--suffix_len&quot; 设置为 11。</p></li></ul><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --Assemble --platform illumina --assembler abyss --filter_length <span class="token number">200</span> --ReadsPath Reads/Illumina --reads1 _1.fastq.gz --reads2 _2.fastq.gz --kmmer <span class="token number">81</span> --threads <span class="token number">4</span> --suffix_len <span class="token number">11</span></pre></td></tr><tr><td data-num="2"></td><td><pre></pre></td></tr><tr><td data-num="3"></td><td><pre><span class="token variable">$pgcgap</span> --Assemble --platform illumina --assembler spades --filter_length <span class="token number">200</span> --ReadsPath Reads/Illumina --reads1 _1.fastq.gz --reads2 _2.fastq.gz --threads <span class="token number">4</span> --suffix_len <span class="token number">11</span></pre></td></tr><tr><td data-num="4"></td><td><pre></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$pgcgap</span> --Assemble --platform illumina --assembler auto --filter_length <span class="token number">200</span> --ReadsPath Reads/Illumina --reads1 _1.fastq.gz --reads2 _2.fastq.gz --kmmer <span class="token number">81</span> --threads <span class="token number">4</span> --suffix_len <span class="token number">11</span></pre></td></tr></table></figure><ul><li><p><strong>Oxford reads 组装</strong></p><p>Oxford nanopore 测序仅产生一个 reads 文件，因此只需要设置 &quot;--reads1&quot; 参数，其值为 &quot;.fasta&quot;。 “--genomeSize” 是预估的基因组大小，用户可以到 NCBI 数据库中查看同物种基因组的大小作为参考，此处设置为 &quot;4.8m&quot;。Reads 文件的后缀名为 &quot;.fasta&quot;，其长度为 6，因此将 &quot;--suffix_len&quot; 设置为 6。</p></li></ul><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --Assemble --platform oxford --filter_length <span class="token number">200</span> --ReadsPath Reads/Oxford --reads1 .fasta --genomeSize <span class="token number">4</span>.8m --threads <span class="token number">4</span> --suffix_len <span class="token number">6</span></pre></td></tr></table></figure><ul><li><strong>PacBio reads 组装</strong></li></ul><p>PacBio 同样只产生一个文件 &quot;pacbio.fastq&quot;，参数设置与 Oxford 类似。此处，文件的后缀名为 &quot;.fastq&quot;，其长度为 6，因此 &quot;--suffix_len&quot; 设置为 6。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --Assemble --platform pacbio --filter_length <span class="token number">200</span> --ReadsPath Reads/PacBio --reads1 .fastq --genomeSize <span class="token number">4</span>.8m --threads <span class="token number">4</span> --suffix_len <span class="token number">6</span></pre></td></tr></table></figure></li><li><p><strong>Example 3:</strong> 基因预测及注释。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --Annotate --scafPath Results/Assembles/Scaf/Illumina --Scaf_suffix -8.fa --genus Escherichia --species “Escherichia coli” --codon <span class="token number">11</span> --threads <span class="token number">4</span></pre></td></tr></table></figure></li><li><p><strong>Example 4:</strong> 构建单拷贝核心蛋白进化树与核心 SNPs 进化树。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment"># Construct phylogenetic tree with FastTree (Quick without best fit model testing)</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$pgcgap</span> --CoreTree --CDsPath Results/Annotations/CDs --AAsPath Results/Annotations/AAs --codon <span class="token number">11</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --fasttree</pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Very slow with best fit model testing, traditional bootstrap)</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$pgcgap</span> --CoreTree --CDsPath Results/Annotations/CDs --AAsPath Results/Annotations/AAs --codon <span class="token number">11</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --bsnum <span class="token number">500</span></pre></td></tr><tr><td data-num="6"></td><td><pre></pre></td></tr><tr><td data-num="7"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Slow with best fit model testing, ultrafast bootstrap)</span></pre></td></tr><tr><td data-num="8"></td><td><pre><span class="token variable">$pgcgap</span> --CoreTree --CDsPath Results/Annotations/CDs --AAsPath Results/Annotations/AAs --codon <span class="token number">11</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --fastboot <span class="token number">1000</span></pre></td></tr></table></figure></li><li><p><strong>Example 5:</strong> 仅构建单拷贝核心蛋白进化树。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment"># Construct phylogenetic tree with FastTree (Quick without best fit model testing)</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$pgcgap</span> --CoreTree --CDsPath NO --AAsPath Results/Annotations/AAs --codon <span class="token number">11</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --fasttree</pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Very slow with best fit model testing, traditional bootstrap)</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$pgcgap</span> --CoreTree --CDsPath NO --AAsPath Results/Annotations/AAs --codon <span class="token number">11</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --bsnum <span class="token number">500</span></pre></td></tr><tr><td data-num="6"></td><td><pre></pre></td></tr><tr><td data-num="7"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Slow with best fit model testing, ultrafast bootstrap)</span></pre></td></tr><tr><td data-num="8"></td><td><pre><span class="token variable">$pgcgap</span> --CoreTree --CDsPath NO --AAsPath Results/Annotations/AAs --codon <span class="token number">11</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --fastboot <span class="token number">1000</span></pre></td></tr></table></figure></li><li><p><strong>Example 6:</strong> 进行泛基因组分析并构建单拷贝核心蛋白进化树。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment"># Construct phylogenetic tree with FastTree (Quick without best fit model testing)</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$pgcgap</span> --Pan --codon <span class="token number">11</span> --identi <span class="token number">95</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --GffPath Results/Annotations/GFF --PanTree --fasttree</pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Very slow with best fit model testing, traditional bootstrap)</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$pgcgap</span> --Pan --codon <span class="token number">11</span> --identi <span class="token number">95</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --GffPath Results/Annotations/GFF --PanTree --bsnum <span class="token number">500</span></pre></td></tr><tr><td data-num="6"></td><td><pre></pre></td></tr><tr><td data-num="7"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Slow with best fit model testing, ultrafast bootstrap)</span></pre></td></tr><tr><td data-num="8"></td><td><pre><span class="token variable">$pgcgap</span> --Pan --codon <span class="token number">11</span> --identi <span class="token number">95</span> --strain_num <span class="token number">6</span> --threads <span class="token number">4</span> --GffPath Results/Annotations/GFF --PanTree --fastboot <span class="token number">1000</span></pre></td></tr></table></figure></li><li><p><strong>Example 7:</strong> 同源蛋白家族聚类分析并构建进化树。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment"># Construct phylogenetic tree with FastTree (Quick without best fit model testing)</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$pgcgap</span> --OrthoF --threads <span class="token number">4</span> --AAsPath Results/Annotations/AAs --fasttree</pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Very slow with best fit model testing, traditional bootstrap)</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$pgcgap</span> --OrthoF --threads <span class="token number">4</span> --AAsPath Results/Annotations/AAs --bsnum <span class="token number">500</span></pre></td></tr><tr><td data-num="6"></td><td><pre></pre></td></tr><tr><td data-num="7"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Slow with best fit model testing, ultrafast bootstrap)</span></pre></td></tr><tr><td data-num="8"></td><td><pre><span class="token variable">$pgcgap</span> --OrthoF --threads <span class="token number">4</span> --AAsPath Results/Annotations/AAs --fastboot <span class="token number">1000</span></pre></td></tr></table></figure></li><li><p><strong>Example 8:</strong> 计算两两基因组之间的平均核苷酸一致性 (ANI)。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --ANI --threads <span class="token number">4</span> --queryL scaf.list --refL scaf.list --ANIO Results/ANI/ANIs --Scaf_suffix .fa</pre></td></tr></table></figure></li><li><p><strong>Example 9:</strong> 通过 MinHash 计算基因组及宏基因组的相似性。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --MASH --scafPath <span class="token operator">&lt;</span><span class="token environment constant">PATH</span><span class="token operator">></span> --Scaf_suffix <span class="token operator">&lt;</span>STRING<span class="token operator">></span></pre></td></tr></table></figure></li><li><p><strong>Example 10:</strong> 对所有基因组进行 COG 注释。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --pCOG --threads <span class="token number">4</span> --strain_num <span class="token number">6</span> --AAsPath Results/Annotations/AAs</pre></td></tr></table></figure></li><li><p><strong>Example 11:</strong> 变异 (SNPs, Indels) 检测与注释，并构建基于参考基因组的 SNPs 进化树。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Very slow with best fit model testing, traditional bootstrap)</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$pgcgap</span> --VAR --threads <span class="token number">4</span> --refgbk /mnt/h/PGCGAP_Examples/Reads/MG1655.gbff --ReadsPath Reads/Illumina --reads1 _1.fastq.gz --reads2 _2.fastq.gz --suffix_len <span class="token number">11</span> --strain_num <span class="token number">6</span> --qualtype sanger --bsnum <span class="token number">500</span></pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Slow with best fit model testing, ultrafast bootstrap)</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$pgcgap</span> --VAR --threads <span class="token number">4</span> --refgbk /mnt/h/PGCGAP_Examples/Reads/MG1655.gbff --ReadsPath Reads/Illumina --reads1 _1.fastq.gz --reads2 _2.fastq.gz --suffix_len <span class="token number">11</span> --strain_num <span class="token number">6</span> --qualtype sanger --fastboot <span class="token number">1000</span></pre></td></tr></table></figure></li><li><p><strong>Example 12:</strong> 从基因组中寻找抗生素抗性基因或毒力基因。</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --AntiRes --scafPath Results/Assembles/Scaf/Illumina --Scaf_suffix -8.fa --threads <span class="token number">6</span> --db ncbi --identity <span class="token number">75</span> --coverage <span class="token number">50</span></pre></td></tr></table></figure></li><li><p><strong>Example 13:</strong> Filter short sequences in the genome and assess the status of the genome</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$pgcgap</span> --ACC --Assess --scafPath Results/Assembles/Scaf/Illumina --Scaf_suffix -8.fa --filter_length <span class="token number">200</span></pre></td></tr></table></figure></li><li><p><strong>Example 14:</strong> Construct a phylogenetic tree based on multiple sequences in one file</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Very slow with best fit model testing, traditional bootstrap)</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$pgcgap</span> --STREE --seqfile proteins.fas --seqtype p --bsnum <span class="token number">500</span> --threads <span class="token number">4</span></pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment"># Construct phylogenetic tree with IQ-TREE (Slow with best fit model testing, ultrafast bootstrap)</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$pgcgap</span> --STREE --seqfile proteins.fas --seqtype p --fastboot <span class="token number">1000</span> --threads <span class="token number">4</span></pre></td></tr></table></figure></li></ul><h2 id="准备输入文件"><a class="anchor" href="#准备输入文件">#</a> 准备输入文件</h2><h3 id="工作目录"><a class="anchor" href="#工作目录">#</a> 工作目录</h3><ul><li>PGCGAP 的运行目录。</li></ul><h3 id="assemble"><a class="anchor" href="#assemble">#</a> Assemble</h3><ul><li>将所有双端 reads、或 PacBio reads 或 Oxford nanopore reads 存于某个目录下 (Default: ./Reads/Illumina/)。</li></ul><h3 id="annotate"><a class="anchor" href="#annotate">#</a> Annotate</h3><ul><li>基因组文件 (完整或不完整) 存放于某个目录下 (Default: Results/Assembles/Scaf/Illumina)。</li></ul><h3 id="ani"><a class="anchor" href="#ani">#</a> ANI</h3><ul><li>QUERY_LIST 和 REFERENCE_LIST 文件，每个文件中含有需要计算的基因组的绝对路径，每个基因组一行 (default: 工作目录下的 &quot;scaf.list&quot;)。若先运行了 “--Assemble” 模块，该文件会自动生成。</li></ul><h3 id="coretree"><a class="anchor" href="#coretree">#</a> CoreTree</h3><ul><li>将所有菌株的氨基酸文件 (后缀名必须为 “.faa”) 和核苷酸文件 (需以 “.ffn” 为后缀) 分别存放于两个目录中 (default: “./Results/Annotations/AAs/” 和 “./Results/Annotations/CDs/”)。“.faa” 和 “.ffn” 文件需要有相同的前缀名字，且 protein IDs 和 gene IDs 需以菌株名开头。建议用 “Prokka” 软件获取输入文件，若已经运行了 “--Annotate” 模块，则该模块的输入文件会自动生产。若 “--CDsPath” 设置为 “NO”，则不需要提供核苷酸文件，但也不会生产核心 SNPs 进化树。</li></ul><h3 id="mash"><a class="anchor" href="#mash">#</a> MASH</h3><ul><li>基因组文件 (完整或不完整) 存放于某个目录下 (Default: Results/Assembles/Scaf/Illumina)。</li></ul><h3 id="orthof"><a class="anchor" href="#orthof">#</a> OrthoF</h3><ul><li>所有菌株的 fasta 格式氨基酸文件 (每个菌株一个文件) 存放于一个目录中 (default: “./Results/Annotations/AAs/”)。 若先运行了 “--Annotate” 模块，该文件会自动生成。</li></ul><h3 id="pan"><a class="anchor" href="#pan">#</a> Pan</h3><ul><li>包含所有菌株 GFF3 文件 (With “.gff” as the suffix) 的目录路径 (default: ./Results/Annotations/GFF/);</li><li>若先运行了 “--Annotate” 模块，上述文件会自动生成。</li></ul><h3 id="pcog"><a class="anchor" href="#pcog">#</a> pCOG</h3><ul><li>存放所有菌株的 fasta 格式氨基酸序列文件 (With “.faa” as the suffix) 的目录路径 (default: ./Results/Annotations/AAs/)。 若先运行了 “--Annotate” 模块，该文件会自动生成。</li></ul><h3 id="var"><a class="anchor" href="#var">#</a> VAR</h3><ul><li>包含所有菌株的 Pair-end reads 的目录路径 (default: ./Reads/Over/under the working directory)。</li><li>fasta 格式或 GenBank 格式的参考基因组的绝对路径 (<strong>必需提供</strong>)。</li></ul><h3 id="antires"><a class="anchor" href="#antires">#</a> AntiRes</h3><ul><li>存放基因组 (complete or draft) 的目录 (Default: Results/Assembles/Scaf/Illumina under the working directory).</li></ul><h3 id="stree"><a class="anchor" href="#stree">#</a> STREE</h3><ul><li>Multiple-FASTA sequences in a file, can be Protein, DNA and Codons.</li></ul><h2 id="输出文件解读"><a class="anchor" href="#输出文件解读">#</a> 输出文件解读</h2><h3 id="assemble-2"><a class="anchor" href="#assemble-2">#</a> Assemble</h3><ul><li><p><strong>Results/Assembles/Illumina/</strong><br><br>Directories contain Illumina assembly files and information of each strain.</p></li><li><p><strong>Results/Assembles/PacBio/</strong><br><br>Directories contain PacBio assembly files and information of each strain.</p></li><li><p><strong>Results/Assembles/Oxford/</strong><br><br>Directories contain ONT assembly files and information of each strain.</p></li><li><p><strong>Results/Assembles/Hybrid/</strong><br><br>Directory contains hybrid assembly files of the short reads and long reads of the same strain.</p></li><li><p><strong>Results/Assembles/Scaf/Illumina</strong><br><br>Directory contains Illumina contigs/scaffolds of all strains. &quot;*.filtered.fas&quot; is the genome after excluding short sequences. &quot;*.prefilter.stats&quot; describes the stats of the genome before filtering, and &quot;*.filtered.stats&quot; describes the stats of the genome after filtering.</p></li><li><p><strong>Results/Assembles/Scaf/Oxford</strong><br><br>Directory contains Oxford nanopore contigs/scaffolds of all strains.</p></li><li><p><strong>Results/Assembles/Scaf/PacBio</strong><br><br>Directory contains PacBio contigs/scaffolds of all strains.</p></li></ul><h3 id="annotate-2"><a class="anchor" href="#annotate-2">#</a> Annotate</h3><ul><li><p><strong>Results/Annotations/*_annotation</strong><br><br>directories contain <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL3Byb2trYT9fYmxhbms=">annotation files</span> of each strain.</p></li><li><p><strong>Results/Annotations/AAs</strong><br><br>Directory contain amino acids sequences of all strains.</p></li><li><p><strong>Results/Annotations/CDs</strong><br><br>Directory contain nucleotide sequences of all strains.</p></li><li><p><strong>Results/Annotations/GFF</strong><br><br>Directory contain the master annotation of all strains in GFF3 format.</p></li></ul><h3 id="ani-2"><a class="anchor" href="#ani-2">#</a> ANI</h3><ul><li><p><strong>Results/ANI/ANIs</strong><br><br>The file contains comparation information of genome pairs. The document is composed of five columns, each of which represents query genome, reference genome, ANI value, count of bidirectional fragment mappings, total query fragments.</p></li><li><p><strong>Results/ANI/ANIs.matrix</strong><br><br>file with identity values arranged in a <span class="exturl" data-url="aHR0cHM6Ly93d3cubW90aHVyLm9yZy93aWtpL1BoeWxpcC1mb3JtYXR0ZWRfZGlzdGFuY2VfbWF0cml4P19ibGFuaw==">phylip-formatted lower triangular matrix</span>.</p></li><li><p><strong>Results/ANI/ANIs.heatmap</strong><br><br>An ANI matrix of all strains.</p></li><li><p><strong>Results/ANI/ANI_matrix.pdf</strong><br><br>The heatmap plot of &quot;ANIs.heatmap&quot;.</p></li></ul><h3 id="mash-2"><a class="anchor" href="#mash-2">#</a> MASH</h3><ul><li><p><strong>Results/MASH/MASH</strong><br><br>The pairwise distance between pair genomes, each column represents Reference-ID, Query-ID, Mash-distance, P-value, and Matching-hashes, respectively.</p></li><li><p><strong>Results/MASH/MASH2</strong><br><br>The pairwise similarity between pair genomes, each column represents Reference-ID, Query-ID, similarity, P-value, and Matching-hashes, respectively.</p></li><li><p><strong>Results/MASH/MASH.heatmap</strong><br><br>A similarity matrix of all genomes.</p></li><li><p><strong>Results/MASH/MASH_matrix.pdf</strong><br><br>A heat map plot of &quot;MASH.heatmap&quot;.</p></li></ul><h3 id="coretree-2"><a class="anchor" href="#coretree-2">#</a> CoreTree</h3><ul><li><p><strong>Results/CoreTrees/ALL.core.protein.fasta</strong><br><br>Concatenated and aligned sequences file of single-copy core proteins.</p></li><li><p><strong>Results/CoreTrees/ALL.core.protein.nwk</strong><br><br>The phylogenetic tree file of single-copy core proteins for all strains constructed by FastTree.</p></li><li><p><strong>Results/CoreTrees/ALL.core.protein.fasta.gb.treefile</strong><br><br>The phylogenetic tree file of single-copy core proteins for all strains constructed by IQ-TREE.</p></li><li><p><strong>Results/CoreTrees/faa2ffn/ALL.core.nucl.fasta</strong><br><br>Concatenated and aligned sequences file of single-copy core genes.</p></li><li><p><strong>Results/CoreTrees/ALL.core.snp.fasta</strong><br><br>Core SNPs of single-copy core genes in fasta format.</p></li><li><p><strong>Results/CoreTrees/ALL.core.snp.nwk</strong><br><br>The phylogenetic tree file of SNPs of single-copy core genes for all strains constructed by FastTree.</p></li><li><p><strong>Results/CoreTrees/ALL.core.snp.fasta.gb.treefile</strong><br><br>The phylogenetic tree file of SNPs of single-copy core genes for all strains constructed by IQ-TREE</p></li><li><p><strong>Results/CoreTrees/&quot;Other_files&quot;</strong><br><br>Intermediate directories and files.</p></li></ul><h3 id="orthof-2"><a class="anchor" href="#orthof-2">#</a> OrthoF</h3><ul><li><p><strong>Results/OrthoFinder/Results_orthoF</strong><br><br>Same as <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL2RhdmlkZW1tcy9PcnRob0ZpbmRlcj9fYmxhbms=">OrthoFinder</span> outputs.</p></li><li><p><strong>Results/OrthoFinder/Results_orthoF/Single_Copy_Orthologue_Tree/</strong><br><br>Directory contains Phylogenetic tree files based on Single Copy Orthologue sequences.<br></p></li><li><p><strong>Results/OrthoFinder/Results_orthoF/Single_Copy_Orthologue_Tree/Single.Copy.Orthologue.nwk</strong><br><br>Phylogenetic tree constructed by FastTree.<br></p></li><li><p><strong>Results/OrthoFinder/Results_orthoF/Single_Copy_Orthologue_Tree/Single.Copy.Orthologue.fasta.gb.treefile</strong><br><br>Phylogenetic tree constructed by IQ-TREE.<br></p></li></ul><h3 id="pan-2"><a class="anchor" href="#pan-2">#</a> Pan</h3><ul><li><p><strong>Results/PanGenome/Pangenome_Pie.pdf</strong><br><br>A 3D pie chart and a fan chart of the breakdown of genes and the number of isolates they are present in.</p></li><li><p><strong>Results/PanGenome/pangenome_frequency.pdf</strong><br><br>A graph with the frequency of genes versus the number of genomes.</p></li><li><p><strong>Results/PanGenome/Pangenome_matrix.pdf</strong><br><br>A figure showing the tree compared to a matrix with the presence and absence of core and accessory genes.</p></li><li><p><strong>Results/PanGenome/Core/Roary.core.protein.fasta</strong><br><br>Alignments of single-copy core proteins called by roary software.</p></li><li><p><strong>Results/PanGenome/Core/Roary.core.protein.nwk</strong><br><br>A phylogenetic tree of Roary.core.protein.fasta constructed by FastTree.</p></li><li><p><strong>Results/PanGenome/Core/Roary.core.protein.fasta.gb.treefile</strong><br><br>A phylogenetic tree of Roary.core.protein.fasta constructed by IQ-TREE.</p></li><li><p><strong>Results/PanGenome/Other_files</strong><br><br>see <span class="exturl" data-url="aHR0cHM6Ly9zYW5nZXItcGF0aG9nZW5zLmdpdGh1Yi5pby9Sb2FyeS8/X2JsYW5r">roary</span> outputs.</p></li></ul><h3 id="pcog-2"><a class="anchor" href="#pcog-2">#</a> pCOG</h3><ul><li><p><strong>*.COG.xml, *.2gi.table, *.2id.table, *.2Sid.table</strong><br><br>Intermediate files.</p></li><li><p><strong>*.2Scog.table</strong><br><br>The super COG table of each strain.</p></li><li><p><strong>*.2Scog.table.pdf</strong><br><br>A plot of super COG table in pdf format.</p></li><li><p><strong>All_flags_relative_abundances.table</strong><br>A table containing the relative abundance of each flag for all strains.</p></li></ul><h3 id="var-2"><a class="anchor" href="#var-2">#</a> VAR</h3><ul><li><p><strong>Results/Variants/directory-named-in-strains</strong><br><br>directories containing substitutions (snps) and insertions/deletions (indels) of each strain. See <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL3NuaXBweT9fYmxhbms=">Snippy</span> outputs for detail.</p></li><li><p><strong>Results/Variants/Core</strong><br><br>The directory containing SNP phylogeny files.</p><ul><li><strong>core.aln</strong> : A core SNP alignment includes only SNP sites.</li><li><strong>core.full.aln</strong> : A whole genome SNP alignment (includes invariant sites).</li><li><strong>core.*.treefile</strong> : Phylogenetic tree of the core SNP alignment based on the best-fit model of evolution selected using IQ-TREE (ignoring possible recombination).</li><li><strong>gubbins.core.full.node_labelled.final_tree.tre</strong> : Phylogenetic tree of the whole genome SNP alignment constructed with <strong>gubbins</strong> (get rid of recombination).</li></ul></li></ul><h3 id="antires-2"><a class="anchor" href="#antires-2">#</a> AntiRes</h3><ul><li><strong>Results/AntiRes/*.tab</strong> : Screening results of each strain.</li><li><strong>Results/AntiRes/summary.txt</strong> : A matrix of gene presence/absence for all strains.</li></ul><h3 id="stree-2"><a class="anchor" href="#stree-2">#</a> STREE</h3><ul><li><strong>Results/STREE/*.aln</strong> : Aligned sequences.</li><li><strong>Results/STREE/*.aln.gb</strong> : Conserved blocks of the sequences.</li><li><strong>Results/STREE/*.aln.gb.treefile</strong> : The final phylogenetic tree.</li></ul><h2 id="使用许可"><a class="anchor" href="#使用许可">#</a> 使用许可</h2><p>PGCGAP 不可商用，其它情况可使用免费 (licensed under GPLv3)。</p><h2 id="反馈与提问"><a class="anchor" href="#反馈与提问">#</a> 反馈与提问</h2><p>若有问题，可在 <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL2xpYW9jaGVubGFucnVvL3BnY2dhcC9pc3N1ZXM/X2JsYW5r">issues page</span> 提出或通过邮件咨询 <span class="exturl" data-url="bWFpbHRvOmxpYW9jaGVubGFucnVvQHdlYm1haWwuaHphdS5lZHUuY24=">liaochenlanruo@webmail.hzau.edu.cn</span>。</p><h2 id="引用"><a class="anchor" href="#引用">#</a> 引用</h2><ul><li><p>If you use this software please cite: Liu H, Xin B, Zheng J, Zhong H, Yu Y, Peng D, Sun M. Build a bioinformatics analysis platform and apply it to routine analysis of microbial genomics and comparative genomics. <em>Protocol exchange</em>, 2020. DOI: <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjIxMjAzL3JzLjIuMjEyMjQvdjU=">10.21203/rs.2.21224/v5</span></p></li><li><p>If you use &quot;--Assemble&quot;, please also cite one or two of <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL09wZW5HZW5lL2Zhc3Rw">Fastp</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjExMDEvZ3IuMjE0MzQ2LjExNg==">ABySS</span>, <span class="exturl" data-url="aHR0cDovL2xpbmsuc3ByaW5nZXIuY29tL2NoYXB0ZXIvMTAuMTAwNyUyRjk3OC0zLTY0Mi0zNzE5NS0wXzEz">SPAdes</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjExMDEvZ3IuMjE1MDg3LjExNg==">Canu</span>, or <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEzNzEvam91cm5hbC5wY2JpLjEwMDU1OTU=">Unicycler</span>.</p></li><li><p>If you use &quot;--Annotate&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly93d3cucGl4aXYubmV0L21lbWJlcl9pbGx1c3QucGhwP21vZGU9bWVkaXVtJmFtcDtpbGx1c3RfaWQ9MjQ2NDIwNjM=">Prokka</span>.</p></li><li><p>If you use &quot;--CoreTree&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvYmlvaW5mb3JtYXRpY3MvYnRsMTU4">CD-HIT</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbmFyL2drZjQzNg==">MAFFT</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbmFyL2drbDMxNQ==">PAL2NAL</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvYmlvaW5mb3JtYXRpY3MvYnRwMzQ4">trimAL</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEzNzEvam91cm5hbC5wb25lLjAwMDk0OTA=">FastTree</span> or <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbW9sYmV2L21zYWEwMTU=">IQ-TREE</span>, and <span class="exturl" data-url="aHR0cHM6Ly9keC5kb2kub3JnLzEwLjEwOTklMkZtZ2VuLjAuMDAwMDU2">SNP-sites</span>.</p></li><li><p>If you use &quot;--Pan&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly9keC5kb2kub3JnLzEwLjEwOTMlMkZiaW9pbmZvcm1hdGljcyUyRmJ0djQyMQ==">Roary</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbmFyL2drZjQzNg==">MAFFT</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvYmlvaW5mb3JtYXRpY3MvYnRwMzQ4">trimAL</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEzNzEvam91cm5hbC5wb25lLjAwMDk0OTA=">FastTree</span> or <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbW9sYmV2L21zYWEwMTU=">IQ-TREE</span>.</p></li><li><p>If you use &quot;--OrthoF&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly9keC5kb2kub3JnLzEwLjExODYlMkZzMTMwNTktMDE5LTE4MzIteQ==">OrthoFinder</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvYmlvaW5mb3JtYXRpY3MvYnRwMzQ4">trimAL</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEzNzEvam91cm5hbC5wb25lLjAwMDk0OTA=">FastTree</span> or <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbW9sYmV2L21zYWEwMTU=">IQ-TREE</span>.</p></li><li><p>If you use &quot;--ANI&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly9keC5kb2kub3JnLzEwLjEwMzglMkZzNDE0NjctMDE4LTA3NjQxLTk=">fastANI</span>.</p></li><li><p>If you use &quot;--MASH&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly9keC5kb2kub3JnLzEwLjExODYlMkZzMTMwNTktMDE2LTA5OTcteA==">Mash</span>.</p></li><li><p>If you use &quot;--VAR&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL25ham9zaGkvc2lja2xl">Sickle</span>, <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL3NuaXBweQ==">Snippy</span>, <span class="exturl" data-url="aHR0cHM6Ly9keC5kb2kub3JnLzEwLjEwOTMlMkZuYXIlMkZna3UxMTk2">Gubbins</span>, <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbW9sYmV2L21zYWEwMTU=">IQ-TREE</span>, and <span class="exturl" data-url="aHR0cHM6Ly9keC5kb2kub3JnLzEwLjQxNjElMkZmbHkuMTk2OTU=">SnpEff</span>.</p></li><li><p>If you use &quot;--AntiRes&quot;, please also cite <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL2FicmljYXRl">Abricate</span> and the corresponding database you used: <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wbWMvYXJ0aWNsZXMvUE1DNjgxMTQxMA==">NCBI AMRFinderPlus</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjc3ODk3MDU=">CARD</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjI3ODI0ODc=">Resfinder</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjQxNDU1MzI=">ARG-ANNOT</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjY1Nzg1NTk=">VFDB</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wdWJtZWQvMjQ3NzcwOTI=">PlasmidFinder</span>, <span class="exturl" data-url="aHR0cHM6Ly93d3cubmNiaS5ubG0ubmloLmdvdi9wbWMvYXJ0aWNsZXMvUE1DNTM0MzEzNi8=">EcOH</span>, or <span class="exturl" data-url="aHR0cHM6Ly9hY2FkZW1pYy5vdXAuY29tL25hci9hcnRpY2xlLzQ4L0QxL0Q1NjEvNTYyNDk3Mw==">MEGARES 2.00</span>.</p></li><li><p>If you use &quot;--STREE&quot;, please also cite <span class="exturl" data-url="aHR0cDovL2V1cm9wZXBtYy5vcmcvYWJzdHJhY3QvTUVELzMwOTc2Nzkz">Muscle</span> , <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvYmlvaW5mb3JtYXRpY3MvYnRwMzQ4">trimAL</span>, and <span class="exturl" data-url="aHR0cHM6Ly9kb2kub3JnLzEwLjEwOTMvbW9sYmV2L21zYWEwMTU=">IQ-TREE</span>.</p></li></ul><h2 id="常见问题"><a class="anchor" href="#常见问题">#</a> 常见问题</h2><h3 id="q1-var-founction-ran-failed-to-get-annotated-vcfs-and-core-results"><a class="anchor" href="#q1-var-founction-ran-failed-to-get-annotated-vcfs-and-core-results">#</a> Q1 VAR founction ran failed to get annotated VCFs and Core results</h3><p>Check the log file named in &quot;strain_name.log&quot; under Results/Variants/&lt;strain_name&gt;/ directory. If you find a sentence like &quot;WARNING: All frames are zero! This seems rather odd, please check that 'frame' information in your 'genes' file is accurate.&quot; This is an snpEff error. Users can install JDK8 to solve this problem.</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$conda</span> <span class="token function">install</span> java-jdk<span class="token operator">=</span><span class="token number">8.0</span>.112</pre></td></tr></table></figure><p>Click <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3RzZWVtYW5uL3NuaXBweS9pc3N1ZXMvMjU5P19ibGFuaw==">here</span> for more solutions.</p><h3 id="q2-could-not-determine-version-of-minced-please-install-version-2-or-higher"><a class="anchor" href="#q2-could-not-determine-version-of-minced-please-install-version-2-or-higher">#</a> Q2 Could not determine version of minced please install version 2 or higher</h3><p>When running prokka of Assemble founction, this error could happened, the error message shows as following:</p><figure class="highlight text"><figcaption data-lang="text"></figcaption><table><tr><td data-num="1"></td><td><pre>Error: A JNI error has occurred, please check your installation and try again</pre></td></tr><tr><td data-num="2"></td><td><pre>Exception in thread "main" java.lang.UnsupportedClassVersionError: minced has been compiled by a more recent version of the Java Runtime (class file version 55.0), this version of the Java Runtime only recognizes class file versions up to 52.0</pre></td></tr><tr><td data-num="3"></td><td><pre>	at java.lang.ClassLoader.defineClass1(Native Method)</pre></td></tr><tr><td data-num="4"></td><td><pre>	at java.lang.ClassLoader.defineClass(ClassLoader.java:763)</pre></td></tr><tr><td data-num="5"></td><td><pre>	at java.security.SecureClassLoader.defineClass(SecureClassLoader.java:142)</pre></td></tr><tr><td data-num="6"></td><td><pre>	at java.net.URLClassLoader.defineClass(URLClassLoader.java:468)</pre></td></tr><tr><td data-num="7"></td><td><pre>	at java.net.URLClassLoader.access$100(URLClassLoader.java:74)</pre></td></tr><tr><td data-num="8"></td><td><pre>	at java.net.URLClassLoader$1.run(URLClassLoader.java:369)</pre></td></tr><tr><td data-num="9"></td><td><pre>	at java.net.URLClassLoader$1.run(URLClassLoader.java:363)</pre></td></tr><tr><td data-num="10"></td><td><pre>	at java.security.AccessController.doPrivileged(Native Method)</pre></td></tr><tr><td data-num="11"></td><td><pre>	at java.net.URLClassLoader.findClass(URLClassLoader.java:362)</pre></td></tr><tr><td data-num="12"></td><td><pre>	at java.lang.ClassLoader.loadClass(ClassLoader.java:424)</pre></td></tr><tr><td data-num="13"></td><td><pre>	at sun.misc.Launcher$AppClassLoader.loadClass(Launcher.java:349)</pre></td></tr><tr><td data-num="14"></td><td><pre>	at java.lang.ClassLoader.loadClass(ClassLoader.java:357)</pre></td></tr><tr><td data-num="15"></td><td><pre>	at sun.launcher.LauncherHelper.checkAndLoadMain(LauncherHelper.java:495)</pre></td></tr><tr><td data-num="16"></td><td><pre>[01:09:40] Could not determine version of minced - please install version 2.0 or higher</pre></td></tr></table></figure><p>Users can downgrade the minced to version 0.3 to solve this problem.</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token variable">$conda</span> <span class="token function">install</span> <span class="token assign-left variable">minced</span><span class="token operator">=</span><span class="token number">0.3</span></pre></td></tr></table></figure><p>Click <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL2Jpb2NvbmRhL2Jpb2NvbmRhLXJlY2lwZXMvcHVsbC8xNTQwNz9fYmxhbms=">here</span> for detail informations.</p><h3 id="q3-dyld-library-not-loaded-rpathlibcrypto100dylib"><a class="anchor" href="#q3-dyld-library-not-loaded-rpathlibcrypto100dylib">#</a> Q3 dyld: Library not loaded: @rpath/libcrypto.1.0.0.dylib</h3><p>This error may happen when running function &quot;VAR&quot; on macOS. It is an error of openssl. Users can solve this problem as following:</p><figure class="highlight bash"><figcaption data-lang="bash"></figcaption><table><tr><td data-num="1"></td><td><pre><span class="token comment">#Firstly, install brew if have not installed before</span></pre></td></tr><tr><td data-num="2"></td><td><pre><span class="token variable">$ruby</span> -e <span class="token string">"<span class="token variable"><span class="token variable">$(</span><span class="token function">curl</span> -fsSL https://raw.githubusercontent.com/Homebrew/install/master/install<span class="token variable">)</span></span>"</span></pre></td></tr><tr><td data-num="3"></td><td><pre></pre></td></tr><tr><td data-num="4"></td><td><pre><span class="token comment">#Install openssl with brew</span></pre></td></tr><tr><td data-num="5"></td><td><pre><span class="token variable">$brew</span> <span class="token function">install</span> openssl</pre></td></tr><tr><td data-num="6"></td><td><pre></pre></td></tr><tr><td data-num="7"></td><td><pre><span class="token comment">#Create the soft link for libraries</span></pre></td></tr><tr><td data-num="8"></td><td><pre><span class="token variable">$ln</span> -s /usr/local/opt/openssl/lib/libcrypto.1.0.0.dylib /usr/local/lib/</pre></td></tr><tr><td data-num="9"></td><td><pre></pre></td></tr><tr><td data-num="10"></td><td><pre><span class="token variable">$ln</span> -s /usr/local/opt/openssl/lib/libssl.1.0.0.dylib /usr/local/lib/</pre></td></tr></table></figure><p>Click <span class="exturl" data-url="aHR0cHM6Ly9naXN0LmdpdGh1Yi5jb20vYWtsYXAvZTg4NTcyMWVmMTVjODY2OGVkMGExZGQ2NGQyZWExYTc=">here</span> for more informations</p><h3 id="q4-use-of-uninitialized-value-in-require-at-encodepm-line-61"><a class="anchor" href="#q4-use-of-uninitialized-value-in-require-at-encodepm-line-61">#</a> Q4 Use of uninitialized value in require at <span class="exturl" data-url="aHR0cDovL0VuY29kZS5wbQ==">Encode.pm</span> line 61</h3><p>This warning may happen when running function &quot;Pan&quot;. It is a warning of Roary software.<br>The content of line 61 is &quot;require Encode::ConfigLocal;&quot;. Users can ignore the warning.<br>Click <span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL3Nhbmdlci1wYXRob2dlbnMvUm9hcnkvaXNzdWVzLzMyMw==">here</span> for details.</p><h2 id="updates"><a class="anchor" href="#updates">#</a> Updates</h2><ul><li><p>V1.0.3</p><ul><li>Updated ANI fuction.</li></ul></li><li><p>V1.0.4</p><ul><li>Add parallel for function &quot;COG&quot;.</li><li>Optimized drawing of ANI heat map.</li></ul></li><li><p>V1.0.5</p><ul><li>Bug repair for input of gubbins.</li></ul></li><li><p>V1.0.6</p><ul><li>Modified CoreTree to split protein and SNPs tree constructing.</li></ul></li><li><p>V1.0.7</p><ul><li>Split Assemble and Annotate into two functions.</li><li>Added third generation genome assembly function.</li><li>Changed the default parameters of the CoreTree function (aS 0.8 to 0.7 and aL 0.8 to 0.5).</li><li>Changed the name of function &quot;COG&quot; to &quot;pCOG&quot;.</li><li>Fixed the sorting bug for ANI heat map.</li></ul></li><li><p>V1.0.8</p><ul><li>Add the &quot;MASH&quot; function to compute genome distance and similarity using MinHash.</li></ul></li><li><p>V1.0.9</p><ul><li>The function of constructing a single-copy core protein phylogenetic tree was added to &quot;Pan&quot;.</li><li>Fixed a bug of plot_3Dpie.R, Optimized image display, and a fan chart has been added.</li><li>Fixed a bug for ploting the ANI matrix.</li></ul></li><li><p>V1.0.10</p><ul><li>Add the &quot;AntiRes&quot; function to screening of contigs for antimicrobial and virulence genes.</li></ul></li><li><p>V1.0.11</p><ul><li>Users now can choose &quot;abyss&quot; or &quot;spades&quot; for illumina reads aseembly.</li><li>New support for hybrid assembly of paired-end short reads and long reads.</li><li>Add the selecting of best-fit model of evolution for DNA and protein alignments before constructing a phylogenetic tree.</li><li>Optimized display of help information. Users can check parameters for each modulewith command &quot;pgcgap [Assemble|Annotate|ANI|AntiRes|CoreTree|MASH|OrthoF|Pan|pCOG|VAR]&quot;, and can look up the examples of each module with command &quot;pgcgap Examples&quot;.</li></ul></li><li><p>V1.0.12</p><ul><li>Added automatic mode for illumina genome assembly. First, PGCGAP calls &quot;ABySS&quot; for genome assembly. When the assembled N50 is less than 50,000, it automatically calls &quot;SPAdes&quot; to try multiple parameters for assembly.</li><li>Added ability to filter short sequences of assembled genomes.</li><li>Added function of genome assembly status assessment.</li><li>Modified the drawing script of ANI and MASH modules so that it can automatically adjust the font size according to the number of samples.</li></ul></li><li><p>V1.0.13</p><ul><li>Fixed the &quot;running error&quot; bug of function &quot;Assess&quot; in module &quot;ACC&quot;.</li><li>Added module &quot;STREE&quot; for constructing a phylogenetic tree based on multiple sequences in one file.</li></ul></li><li><p>V1.0.14</p><ul><li>The relative_abundances of flags among strains will not be called while the strain number is less than two.</li></ul></li><li><p>V1.0.15</p><ul><li>When the number of threads set by the user exceeds the number of threads owned by the system, PGCGAP will automatically adjust the number of threads to avoid program crash.</li><li>Add FASTQ preprocessor before Illunima genome assembly: adapter trimming, polyG tail trimming of Illumina NextSeq/NovaSeq reads, quality filtering (Q value filtering, N base filtering, sliding window filtering), length filtering.</li></ul></li><li><p>V1.0.16</p><ul><li>Reduced the number of Racon polishing rounds for better speed performance when peforming genome assembly.</li><li>Force overwriting existing output folder when running &quot;Annotate&quot; analysis to avoid program crash.</li></ul></li><li><p>V1.0.17</p><ul><li>Fixed a bug that the program can not go back to the working directory after genome annotation.</li><li>Added scripts to check if there were single-copy core proteins found while running module &quot;CoreTree&quot;.</li><li>Modified the help message.</li></ul></li><li><p>V1.0.18</p><ul><li>Updated the downloading link of COG database.</li><li>Users can choose the number of threads used for running module &quot;STREE&quot;.</li></ul></li><li><p>V1.0.19</p><ul><li>Can resume from break-point when downloading the COG database.</li><li>Fixed a bug that failed to create multi-level directories.</li></ul></li><li><p>V1.0.20</p><ul><li>Fixed a little bug (path error) of module &quot;VAR&quot;.</li><li>Fixed a little bug of module &quot;CoreTree&quot; to avoid the interference of special characters in sequence ID to the program.</li></ul></li><li><p>V1.0.21</p><ul><li>Change the default search program &quot;blast&quot; to &quot;diamond&quot; of module &quot;OrthoF&quot;.</li><li>Fixed a bug of module &quot;pCOG&quot; to output the right figure.</li></ul></li><li><p>V1.0.22</p><ul><li>The drawing function of module &quot;ANI&quot; and module &quot;MASH&quot; has been enhanced, including automatic adjustment of font size and legend according to the size of the picture.</li><li>Fixed a bug of module &quot;ANI&quot;, that is no heatmap will be drawn when there is &quot;NA&quot; in the ANI matrix in the previous versions.</li><li>When the ANI value or genome similarity is greater than 95%, an asterisk (*) will be drawn in the corresponding cell of the heatmap.</li></ul></li><li><p>V1.0.23</p><ul><li>The &quot;--Assess&quot; function of module &quot;ACC&quot; was enhanced to (1) generate a summary file containing the status of all genomes (before and after the short sequence filtering), (2) auto move the low-quality genomes (that is genomes with N50 length less than 50 k) to a directory, and others to another directory.</li></ul></li><li><p>V1.0.24</p><ul><li>Fixed a little bug of module &quot;Pan&quot; to avoid the interference of special characters (&gt;) in sequence ID to the program.</li></ul></li><li><p>V1.0.25</p><ul><li>Gblocks was used to eliminate poorly aligned positions and divergent regions of an alignment of DNA or protein sequences in module &quot;CoreTree&quot; and &quot;Pan&quot;.</li><li>The parameter &quot;--identi&quot; was added into module &quot;Pan&quot; to allow users to set the minimum percentage identity for blastp.</li></ul></li><li><p>V1.0.26</p><ul><li>Adjusted the font size with the variation of genome number and the string length of the genome name when plotting the heat map of module &quot;ANI&quot; and &quot;MASH&quot;.</li><li>Two heat map are provided, one of which with a star (means the similarity of the two genomes is larger than 95%) and another without a star, when performing the &quot;ANI&quot; and &quot;MASH&quot; analysis.</li></ul></li><li><p>V1.0.27</p><ul><li>The Amino Acid files are no longer needed when performing the Pan-genome analysis with module Pan.</li></ul></li><li><p>V1.0.28</p><ul><li>Users can check and install the latest version of PGCGAP by the command &quot;pgcgap --check-update&quot;.</li><li>Update module Assemble to allow polish after the assembly of PacBio and ONT data.</li><li>Update module pCOG to adjust the latest database of <span class="exturl" data-url="aHR0cHM6Ly9mdHAubmNiaS5uaWguZ292L3B1Yi9DT0cvQ09HMjAyMC9kYXRhL0NPRw==">COG 2020</span>.</li><li>Optimized the drawing and color scheme of the module pCOG.</li><li>Fixed the parameter &quot;CoreTree&quot; in the module Pan to avoid program termination caused by the &quot;&gt;&quot; in non-sequence lines.</li></ul></li><li><p>V1.0.29</p><ul><li>Function added to module OrthoF: Phylogenetic tree can be constructed automatically with the Single Copy Orthologue Sequences called by module OrthoF.</li><li>Fixed the &quot;permission denied&quot; error when moving directories on the WSL platform.</li></ul></li><li><p>V1.0.30</p></li><li><p>Replace Gblocks with trimAL to trim MSA (module CoreTree, Pan, STREE, and OrthoF).</p></li><li><p>Replaced Modeltest-ng and Raxml-ng with IQ-TREE (module CoreTree, Pan OrthoF, and VAR).</p></li><li><p>Added the option of using FastTree to build phylogenetic tree (module CoreTree, Pan, and OrthoF).</p></li><li><p>V1.0.31</p><ul><li>The default replicates for bootstrap testing was set to 500.</li><li>Add the method for phylogenetic tree constructing with ultrafast bootstrap of IQ-TREE.</li><li>Prevent the log from being written to the tree file generated by FastTree.</li></ul></li><li><p>V1.0.32</p><ul><li>A more colorful version, try &quot;pgcgap Examples&quot; to have a look.</li><li>Updated module AntiRes: the parameter --db had been modified to add choices of &quot;all&quot; and &quot;megares&quot;.</li><li>A little optimization of module VAR.</li><li>Replaced conda with mamba to update PGCGAP more quickly.</li></ul></li></ul><div class="tags"><a href="/tags/%E6%AF%94%E8%BE%83%E5%9F%BA%E5%9B%A0%E7%BB%84%E5%AD%A6/" rel="tag"><i class="ic i-tag"></i> 比较基因组学</a></div></div><footer><div class="meta"><span class="item"><span class="icon"><i class="ic i-calendar-check"></i> </span><span class="text">Edited on</span> <time title="Modified: 2022-05-31 09:41:33" itemprop="dateModified" datetime="2022-05-31T09:41:33+08:00">2022-05-31</time> </span><span id="post/848f.html" class="item leancloud_visitors" data-flag-title="PGCGAP中文说明" title="Views"><span class="icon"><i class="ic i-eye"></i> </span><span class="text">Views</span> <span class="leancloud-visitors-count"></span> <span class="text">times</span></span></div><div class="reward"><button><i class="ic i-heartbeat"></i> Donate</button><p>Give me a cup of [coffee]~(￣▽￣)~*</p><div id="qr"><div><img data-src="/images/reward-wepays.jpg" alt="Hualin Liu WeChat Pay"><p>WeChat Pay</p></div><div><img data-src="/images/AliPays.jpg" alt="Hualin Liu Alipay"><p>Alipay</p></div></div></div><div id="copyright"><ul><li class="author"><strong>Post author: </strong>liaochenlanruo <i class="ic i-at"><em>@</em></i>了尘兰若的小坑</li><li class="link"><strong>Post link: </strong><a href="https://liaochenlanruo.gitee.io/post/848f.html" title="PGCGAP中文说明">https://liaochenlanruo.gitee.io/post/848f.html</a></li><li class="license"><strong>Copyright Notice: </strong>All articles in this blog are licensed under <span class="exturl" data-url="aHR0cHM6Ly9jcmVhdGl2ZWNvbW1vbnMub3JnL2xpY2Vuc2VzL2J5LW5jLXNhLzQuMC9kZWVkLnpo"><i class="ic i-creative-commons"><em>(CC)</em></i>BY-NC-SA</span> unless stating additionally.</li></ul></div></footer></article></div><div class="post-nav"><div class="item left"><a href="/post/9092.html" itemprop="url" rel="prev" data-background-image="https:&#x2F;&#x2F;fastly.jsdelivr.net&#x2F;gh&#x2F;liaochenlanruo&#x2F;cdn@master&#x2F;img&#x2F;custom&#x2F;bgs&#x2F;thumb_5.webp" title="Perl获取外部命令执行结果的输出"><span class="type">Previous Post</span> <span class="category"><i class="ic i-flag"></i> IT</span><h3>Perl获取外部命令执行结果的输出</h3></a></div><div class="item right"><a href="/post/d163.html" itemprop="url" rel="next" data-background-image="https:&#x2F;&#x2F;fastly.jsdelivr.net&#x2F;gh&#x2F;liaochenlanruo&#x2F;cdn@master&#x2F;img&#x2F;custom&#x2F;bgssmall&#x2F;bioconda.png" title="Bioconda贡献指南"><span class="type">Next Post</span> <span class="category"><i class="ic i-flag"></i> 生物信息</span><h3>Bioconda贡献指南</h3></a></div></div><div class="wrap" id="comments"></div></div><div id="sidebar"><div class="inner"><div class="panels"><div class="inner"><div class="contents panel pjax" data-title="Contents"><ol class="toc"><li class="toc-item toc-level-2"><a class="toc-link" href="#%E7%AE%80%E4%BB%8B"><span class="toc-number">1.</span> <span class="toc-text">简介</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E5%AE%89%E8%A3%85"><span class="toc-number">2.</span> <span class="toc-text">安装</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E4%BE%9D%E8%B5%96%E5%8C%85"><span class="toc-number">3.</span> <span class="toc-text">依赖包</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#pgcgap%E7%94%A8%E6%B3%95"><span class="toc-number">4.</span> <span class="toc-text">PGCGAP 用法</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E7%A4%BA%E4%BE%8B"><span class="toc-number">5.</span> <span class="toc-text">示例</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E5%87%86%E5%A4%87%E8%BE%93%E5%85%A5%E6%96%87%E4%BB%B6"><span class="toc-number">6.</span> <span class="toc-text">准备输入文件</span></a><ol class="toc-child"><li class="toc-item toc-level-3"><a class="toc-link" href="#%E5%B7%A5%E4%BD%9C%E7%9B%AE%E5%BD%95"><span class="toc-number">6.1.</span> <span class="toc-text">工作目录</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#assemble"><span class="toc-number">6.2.</span> <span class="toc-text">Assemble</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#annotate"><span class="toc-number">6.3.</span> <span class="toc-text">Annotate</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#ani"><span class="toc-number">6.4.</span> <span class="toc-text">ANI</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#coretree"><span class="toc-number">6.5.</span> <span class="toc-text">CoreTree</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#mash"><span class="toc-number">6.6.</span> <span class="toc-text">MASH</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#orthof"><span class="toc-number">6.7.</span> <span class="toc-text">OrthoF</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#pan"><span class="toc-number">6.8.</span> <span class="toc-text">Pan</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#pcog"><span class="toc-number">6.9.</span> <span class="toc-text">pCOG</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#var"><span class="toc-number">6.10.</span> <span class="toc-text">VAR</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#antires"><span class="toc-number">6.11.</span> <span class="toc-text">AntiRes</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#stree"><span class="toc-number">6.12.</span> <span class="toc-text">STREE</span></a></li></ol></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E8%BE%93%E5%87%BA%E6%96%87%E4%BB%B6%E8%A7%A3%E8%AF%BB"><span class="toc-number">7.</span> <span class="toc-text">输出文件解读</span></a><ol class="toc-child"><li class="toc-item toc-level-3"><a class="toc-link" href="#assemble-2"><span class="toc-number">7.1.</span> <span class="toc-text">Assemble</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#annotate-2"><span class="toc-number">7.2.</span> <span class="toc-text">Annotate</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#ani-2"><span class="toc-number">7.3.</span> <span class="toc-text">ANI</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#mash-2"><span class="toc-number">7.4.</span> <span class="toc-text">MASH</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#coretree-2"><span class="toc-number">7.5.</span> <span class="toc-text">CoreTree</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#orthof-2"><span class="toc-number">7.6.</span> <span class="toc-text">OrthoF</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#pan-2"><span class="toc-number">7.7.</span> <span class="toc-text">Pan</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#pcog-2"><span class="toc-number">7.8.</span> <span class="toc-text">pCOG</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#var-2"><span class="toc-number">7.9.</span> <span class="toc-text">VAR</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#antires-2"><span class="toc-number">7.10.</span> <span class="toc-text">AntiRes</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#stree-2"><span class="toc-number">7.11.</span> <span class="toc-text">STREE</span></a></li></ol></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E4%BD%BF%E7%94%A8%E8%AE%B8%E5%8F%AF"><span class="toc-number">8.</span> <span class="toc-text">使用许可</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E5%8F%8D%E9%A6%88%E4%B8%8E%E6%8F%90%E9%97%AE"><span class="toc-number">9.</span> <span class="toc-text">反馈与提问</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E5%BC%95%E7%94%A8"><span class="toc-number">10.</span> <span class="toc-text">引用</span></a></li><li class="toc-item toc-level-2"><a class="toc-link" href="#%E5%B8%B8%E8%A7%81%E9%97%AE%E9%A2%98"><span class="toc-number">11.</span> <span class="toc-text">常见问题</span></a><ol class="toc-child"><li class="toc-item toc-level-3"><a class="toc-link" href="#q1-var-founction-ran-failed-to-get-annotated-vcfs-and-core-results"><span class="toc-number">11.1.</span> <span class="toc-text">Q1 VAR founction ran failed to get annotated VCFs and Core results</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#q2-could-not-determine-version-of-minced-please-install-version-2-or-higher"><span class="toc-number">11.2.</span> <span class="toc-text">Q2 Could not determine version of minced please install version 2 or higher</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#q3-dyld-library-not-loaded-rpathlibcrypto100dylib"><span class="toc-number">11.3.</span> <span class="toc-text">Q3 dyld: Library not loaded: @rpath&#x2F;libcrypto.1.0.0.dylib</span></a></li><li class="toc-item toc-level-3"><a class="toc-link" href="#q4-use-of-uninitialized-value-in-require-at-encodepm-line-61"><span class="toc-number">11.4.</span> <span class="toc-text">Q4 Use of uninitialized value in require at Encode.pm line 61</span></a></li></ol></li><li class="toc-item toc-level-2"><a class="toc-link" href="#updates"><span class="toc-number">12.</span> <span class="toc-text">Updates</span></a></li></ol></div><div class="related panel pjax" data-title="Related"><ul><li><a href="/post/19824.html" rel="bookmark" title="生物信息学1：VMware虚拟机及Bio-linux安装与配置">生物信息学1：VMware虚拟机及Bio-linux安装与配置</a></li><li><a href="/post/9.html" rel="bookmark" title="生物信息学2：VirtualBox虚拟机及Bio-Linux安装">生物信息学2：VirtualBox虚拟机及Bio-Linux安装</a></li><li><a href="/post/10877.html" rel="bookmark" title="生物信息学3：微生物基因组学常用软件安装">生物信息学3：微生物基因组学常用软件安装</a></li><li><a href="/post/30650.html" rel="bookmark" title="根据基因组预测表型 —— traitar的安装与使用">根据基因组预测表型 —— traitar的安装与使用</a></li><li><a href="/post/44606.html" rel="bookmark" title="kSNP3寻找SNPs并构建进化树">kSNP3寻找SNPs并构建进化树</a></li><li><a href="/post/43504.html" rel="bookmark" title="用wget批量下载含有链接的文件/目录">用wget批量下载含有链接的文件/目录</a></li><li><a href="/post/94f2.html" rel="bookmark" title="为PubMed添加功能">为PubMed添加功能</a></li><li><a href="/post/4e1.html" rel="bookmark" title="OrthoFinder2—同源蛋白家族聚类">OrthoFinder2—同源蛋白家族聚类</a></li><li class="active"><a href="/post/848f.html" rel="bookmark" title="PGCGAP中文说明">PGCGAP中文说明</a></li><li><a href="/post/d163.html" rel="bookmark" title="Bioconda贡献指南">Bioconda贡献指南</a></li><li><a href="/post/f6c9.html" rel="bookmark" title="Linux和基因组测序技术初识">Linux和基因组测序技术初识</a></li><li><a href="/post/5908.html" rel="bookmark" title="BtToxin_Digger--A high-throughput Bacillus thuringiensis toxin mining pipeline">BtToxin_Digger--A high-throughput Bacillus thuringiensis toxin mining pipeline</a></li><li><a href="/post/ab1b.html" rel="bookmark" title="宏基因组分析流程及代码">宏基因组分析流程及代码</a></li><li><a href="/post/4cf4.html" rel="bookmark" title="Use microeco分析扩增子数据">Use microeco分析扩增子数据</a></li><li><a href="/post/1755.html" rel="bookmark" title="扩增子分析--计算随机过程和决定性过程比例">扩增子分析--计算随机过程和决定性过程比例</a></li><li><a href="/post/e7e9.html" rel="bookmark" title="根据assession number批量从NCB下载数据">根据assession number批量从NCB下载数据</a></li><li><a href="/post/e922.html" rel="bookmark" title="Swissprot数据库的本地化与序列比对并与其他数据库快速mapping">Swissprot数据库的本地化与序列比对并与其他数据库快速mapping</a></li><li><a href="/post/412a.html" rel="bookmark" title="构建样本vs基因矩阵">构建样本vs基因矩阵</a></li><li><a href="/post/9052.html" rel="bookmark" title="CAZy碳水化合物活性酶预测">CAZy碳水化合物活性酶预测</a></li><li><a href="/post/5324.html" rel="bookmark" title="原核生物基因岛预测">原核生物基因岛预测</a></li><li><a href="/post/d9f9.html" rel="bookmark" title="按照Contig切割GenBank文件">按照Contig切割GenBank文件</a></li><li><a href="/post/3a23.html" rel="bookmark" title="从GenBank文件中提取Features">从GenBank文件中提取Features</a></li><li><a href="/post/fb68.html" rel="bookmark" title="使用EffectiveT3预测微生物中的III型分泌系统效应蛋白">使用EffectiveT3预测微生物中的III型分泌系统效应蛋白</a></li><li><a href="/post/c886.html" rel="bookmark" title="SignalP+TMHMM预测微生物分泌蛋白">SignalP+TMHMM预测微生物分泌蛋白</a></li><li><a href="/post/6a93.html" rel="bookmark" title="在FASTA文件中搜索完全匹配的短序列">在FASTA文件中搜索完全匹配的短序列</a></li><li><a href="/post/ea71.html" rel="bookmark" title="在Linux服务器中安装网页版Rstudio">在Linux服务器中安装网页版Rstudio</a></li><li><a href="/post/a713.html" rel="bookmark" title="R语言安装依赖包错误集锦">R语言安装依赖包错误集锦</a></li><li><a href="/post/d444.html" rel="bookmark" title="NCBI上传基因簇之tbl2asn的使用">NCBI上传基因簇之tbl2asn的使用</a></li><li><a href="/post/18ea.html" rel="bookmark" title="利用NCycDB数据库从宏基因组中预测氮循环基因">利用NCycDB数据库从宏基因组中预测氮循环基因</a></li><li><a href="/post/36cc.html" rel="bookmark" title="R语言绘制分组条形图">R语言绘制分组条形图</a></li><li><a href="/post/1551.html" rel="bookmark" title="宏病毒组分析流程1-VirSorter2">宏病毒组分析流程1-VirSorter2</a></li><li><a href="/post/8fb3.html" rel="bookmark" title="利用GTDB-TK对细菌和古菌基因组进行物种分类">利用GTDB-TK对细菌和古菌基因组进行物种分类</a></li><li><a href="/post/85d7.html" rel="bookmark" title="计算蛋白质等电点并绘制全局pI图">计算蛋白质等电点并绘制全局pI图</a></li><li><a href="/post/92eb.html" rel="bookmark" title="使用DeepARG预测抗生素抗性基因ARGs">使用DeepARG预测抗生素抗性基因ARGs</a></li><li><a href="/post/74ad.html" rel="bookmark" title="Nanopore基因组组装">Nanopore基因组组装</a></li><li><a href="/post/2581.html" rel="bookmark" title="使用eggNOG-mapper注释基因组">使用eggNOG-mapper注释基因组</a></li><li><a href="/post/d862.html" rel="bookmark" title="PGCGAP - the Prokaryotic Genomics and Comparative Genomics Analysis Pipeline v1.0.33">PGCGAP - the Prokaryotic Genomics and Comparative Genomics Analysis Pipeline v1.0.33</a></li><li><a href="/post/e6ae.html" rel="bookmark" title="利用PGCGAP根据ids提取序列信息">利用PGCGAP根据ids提取序列信息</a></li><li><a href="/post/213.html" rel="bookmark" title="多种方法批量下载NCBI基因组">多种方法批量下载NCBI基因组</a></li><li><a href="/post/33b.html" rel="bookmark" title="代谢组相关软件的安装及使用">代谢组相关软件的安装及使用</a></li><li><a href="/post/377e.html" rel="bookmark" title="使用xcms3处理和分析LC-MS数据">使用xcms3处理和分析LC-MS数据</a></li><li><a href="/post/df57.html" rel="bookmark" title="用metid构建代谢组数据库">用metid构建代谢组数据库</a></li><li><a href="/post/0.html" rel="bookmark" title="基于TidyMass的非靶向代谢组学分析">基于TidyMass的非靶向代谢组学分析</a></li><li><a href="/post/29bf.html" rel="bookmark" title="批量下载某研究方向重要文献">批量下载某研究方向重要文献</a></li><li><a href="/post/c617.html" rel="bookmark" title="16S扩增子测序及宏基因组测序">16S扩增子测序及宏基因组测序</a></li><li><a href="/post/1177.html" rel="bookmark" title="用MaAsLin2包做微生物组多变量线性模型关联分析">用MaAsLin2包做微生物组多变量线性模型关联分析</a></li><li><a href="/post/7449.html" rel="bookmark" title="NCBI上传基因簇之table2asn的使用">NCBI上传基因簇之table2asn的使用</a></li></ul></div><div class="overview panel" data-title="Overview"><div class="author" itemprop="author" itemscope itemtype="http://schema.org/Person"><img class="image" itemprop="image" alt="Hualin Liu" data-src="/images/head.jpg"><p class="name" itemprop="name">Hualin Liu</p><div class="description" itemprop="description">分享微生物生物信息学分析方法，欢迎加入QQ群交流945751012，不接受群内广告！</div></div><nav class="state"><div class="item posts"><a href="/archives/"><span class="count">130</span> <span class="name">posts</span></a></div><div class="item categories"><a href="/categories/"><span class="count">12</span> <span class="name">categories</span></a></div><div class="item tags"><a href="/tags/"><span class="count">54</span> <span class="name">tags</span></a></div></nav><div class="social"><span class="exturl item github" data-url="aHR0cDovL2dpdGh1Yi5jb20vbGlhb2NoZW5sYW5ydW8=" title="http:&#x2F;&#x2F;github.com&#x2F;liaochenlanruo"><i class="ic i-github"></i></span> <span class="exturl item zhihu" data-url="aHR0cHM6Ly93d3cuemhpaHUuY29tL3Blb3BsZS9zdWFuLWxhLW1pYW4tcGlhbi10YW5nLWNpdW0=" title="https:&#x2F;&#x2F;www.zhihu.com&#x2F;people&#x2F;suan-la-mian-pian-tang-cium"><i class="ic i-zhihu"></i></span> <span class="exturl item weibo" data-url="aHR0cHM6Ly93ZWliby5jb20vdS8yNzE1MjEwOTIz" title="https:&#x2F;&#x2F;weibo.com&#x2F;u&#x2F;2715210923"><i class="ic i-weibo"></i></span> <span class="exturl item about" data-url="aHR0cHM6Ly93d3cubGlhb2NoZW5sYW5ydW8uZnVuL2Fib3V0" title="https:&#x2F;&#x2F;www.liaochenlanruo.fun&#x2F;about"><i class="ic i-address-card"></i></span> <span class="exturl item email" data-url="bWFpbHRvOmxpYW9jaGVubGFucnVvQHdlYm1haWwuaHphdS5lZHUuY24=" title="mailto:liaochenlanruo@webmail.hzau.edu.cn"><i class="ic i-envelope"></i></span> <span class="exturl item douban" data-url="aHR0cHM6Ly93d3cuZG91YmFuLmNvbS9wZW9wbGUvbGlhb2NoZW5sYW5ydW8vP19pPTczNDYxMjVRalo1QjZWLw==" title="https:&#x2F;&#x2F;www.douban.com&#x2F;people&#x2F;liaochenlanruo&#x2F;?_i&#x3D;7346125QjZ5B6V&#x2F;"><i class="ic i-douban"></i></span> <span class="exturl item google" data-url="aHR0cHM6Ly9zY2hvbGFyLmdvb2dsZS5jby5qcC9jaXRhdGlvbnM/dXNlcj1YYW9STk1RQUFBQUomaGw9amE=" title="https:&#x2F;&#x2F;scholar.google.co.jp&#x2F;citations?user&#x3D;XaoRNMQAAAAJ&amp;hl&#x3D;ja"><i class="ic i-google"></i></span></div><ul class="menu"><li class="item"><a href="/" rel="section"><i class="ic i-home"></i>Home</a></li><li class="item dropdown"><a href="javascript:void(0);"><i class="ic i-user"></i>About</a><ul class="submenu"><li class="item"><a href="/about/" rel="section"><i class="ic i-user"></i>About me</a></li><li class="item"><a href="/comment/" rel="section"><i class="ic i-envelope"></i>Comment</a></li></ul></li><li class="item dropdown"><a href="javascript:void(0);"><i class="ic i-feather"></i>Posts</a><ul class="submenu"><li class="item"><a href="/archives/" rel="section"><i class="ic i-list-alt"></i>Archives</a></li><li class="item"><a href="/categories/" rel="section"><i class="ic i-th"></i>Categories</a></li><li class="item"><a href="/tags/" rel="section"><i class="ic i-tags"></i>Tags</a></li></ul></li><li class="item dropdown"><a href="javascript:void(0);"><i class="ic i-list-ol"></i>List</a><ul class="submenu"><li class="item"><a href="/photos/" rel="section"><i class="ic i-person"></i>Photos</a></li><li class="item"><a href="/skip/gallery/" rel="section"><i class="ic i-cloud"></i>Gallery</a></li><li class="item"><a href="/skip/box/" rel="section"><i class="ic i-magic"></i>Box</a></li><li class="item"><span class="exturl" data-url="aHR0cHM6Ly9saWFvY2hlbmxhbnJ1by5mdW4vcmVmZXJlbmNlLw=="><i class="ic i-th"></i>cheatsheets</span></li></ul></li><li class="item"><a href="/donate/" rel="section"><i class="ic i-coffee"></i>Donate</a></li><li class="item"><a href="/friends/" rel="section"><i class="ic i-heart"></i>Friends</a></li><li class="item"><span class="exturl" data-url="aHR0cHM6Ly9saWFvY2hlbmxhbnJ1by5naXRodWIuaW8="><i class="ic i-external-link-alt"></i>中文版</span></li><li class="item"><span class="exturl" data-url="aHR0cHM6Ly90cmF2ZWxsaW5ncy5saW5r"><i class="ic i-forward"></i>Travelling</span></li></ul></div></div></div><ul id="quick"><li class="prev pjax"><a href="/post/9092.html" rel="prev" title="Previous Post"><i class="ic i-chevron-left"></i></a></li><li class="up"><i class="ic i-arrow-up"></i></li><li class="down"><i class="ic i-arrow-down"></i></li><li class="next pjax"><a href="/post/d163.html" rel="next" title="Next Post"><i class="ic i-chevron-right"></i></a></li><li class="percent"></li></ul></div></div><div class="dimmer"></div></div></main><footer id="footer"><div class="inner"><div class="widgets"></div><div class="status"><div class="copyright">&copy; 2018 – <span itemprop="copyrightYear">2024</span> <span class="with-love"><i class="ic i-sakura rotate"></i> </span><span class="author" itemprop="copyrightHolder">Hualin Liu @ ResearchGo</span></div><div class="count"><span class="post-meta-item-icon"><i class="ic i-chart-area"></i> </span><span title="Symbols count total">788k words</span> <span class="post-meta-divider">|</span> <span class="post-meta-item-icon"><i class="ic i-coffee"></i> </span><span title="Reading time total">11:57</span><script async src="//busuanzi.ibruce.info/busuanzi/2.3/busuanzi.pure.mini.js"></script><span class="post-meta-divider">|</span> <span class="post-meta-item-icon"><i class="ic i-info-circle"></i> </span><span style="padding-left:1px"><span id="busuanzi_container_site_pv">Total visits: <span id="busuanzi_value_site_pv" style="color:#f99"></span> times</span> <span class="post-meta-divider">|</span> <span class="post-meta-item-icon"><i class="ic i-user"></i> </span><span style="padding-left:1px"><span id="busuanzi_container_site_uv">Total visitors: <span id="busuanzi_value_site_uv" style="color:#f99"></span> people</span></span></span></div><div class="powered-by">Powered by <span class="exturl" data-url="aHR0cHM6Ly9oZXhvLmlv">Hexo</span> & Theme.<span class="exturl" data-url="aHR0cHM6Ly9naXRodWIuY29tL2FtZWhpbWUvaGV4by10aGVtZS1zaG9rYQ==">Shoka</span></div></div></div></footer></div><script data-config type="text/javascript">var LOCAL={path:"post/848f.html",favicon:{show:"（●´3｀●）Goooood",hide:"(´Д｀)Booooom"},search:{placeholder:"Search for Posts",empty:"We didn't find any results for the search: ${query}",stats:"${hits} results found in ${time} ms"},valine:!0,fancybox:!0,copyright:'Copied to clipboard successfully! <br> All articles in this blog are licensed under <i class="ic i-creative-commons"></i>BY-NC-SA.',ignores:[function(e){return e.includes("#")},function(e){return new RegExp(LOCAL.path+"$").test(e)}]}</script><script src="https://cdn.polyfill.io/v2/polyfill.js"></script><script src="//cdn.jsdelivr.net/combine/npm/pace-js@1.0.2/pace.min.js,npm/pjax@0.2.8/pjax.min.js,npm/whatwg-fetch@3.4.0/dist/fetch.umd.min.js,npm/animejs@3.2.0/lib/anime.min.js,npm/algoliasearch@4/dist/algoliasearch-lite.umd.js,npm/instantsearch.js@4/dist/instantsearch.production.min.js,npm/lozad@1/dist/lozad.min.js,npm/quicklink@2/dist/quicklink.umd.js"></script><script src="/js/app.js?v=0.2.5"></script><script data-pjax>var _hmt=_hmt||[];!function(){var e=document.createElement("script");e.src="https://hm.baidu.com/hm.js?d211830a796ceefb38a141fde14e4d5a";var t=document.getElementsByTagName("script")[0];t.parentNode.insertBefore(e,t)}()</script></body></html>